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ABSTRACT 


A time-dependent mixing height model capable of being used 
on a minicomputer system is developed. The main source of energy 
considered in the model is the net surface radiation which is 
calculated as a function of time by a radiation model... A relation- 
ship is derived to transform energy to an equivalent area on a 
thermodynamic diagram. This area is analyzed using basic geometry 
to produce a value of the mixing height. 

It is assumed that thermal turbulence driven by the 
surface radiation excess is the predominant process in the develop- 
ment of the mixed layer. It is assumed also that the morning lapse 
rate of temperature is linear and that the vertical temperature 
gradient in the mixed layers remains constant with time, 

The model requires input data from the surface and the 
850 and 700 hectopascal levels, Required surface data include the 
location's latitude, longitude, elevation above Cees sea level as 
well as the surface albedo and a forecast of the minimum and maximum 
surface temperatures for the day. The heights and temperatures of 
the 850 and 700 hectopascal levels at 1200(GMT) and OOQCO(GMT) are 
used as input. Finally, the cloud albedoes are required at three- 
hourly intervals from 1300(GMT) to O300(GMT) the following day. 

The model presented in this study adequately describes 


the variations in height of the mixed layer for the period from 
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sunrise to late afternoon. A comparison with Holzworth's technique 
(1964), based on test samples of minisonde data in Central Alberta, 
demonstrates a marked improvement in the prediction of mixing heights 
by means’of the new model. Because the input data are readily avail- 
able, the new model is very attractive for operational use where a 
good first approximation of the mixing depth is required. This 
method can also be coupled with a pollution concentration model to 
provide pollution potentials over areas where observed data are 


scarce. 
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CHAPTER 1 


INTRODUCTION 


Pe) the Concept of Mixing Depth 


The mixing layer of the atmosphere comprises the lower 
few kilometers Sree cain that region where surface effects 
predominate. As will be seen later, the mixing layer is synonymous 
with what is called the Planetary Boundary Layer (PBL). The actual 
definition of the layer varies from person to person depending on 
which meteorological parameters are being considered. Summers (1965) 
considered the mixing depth to be the distance from the ground to 
the top /of \the apse layer which formed over an urban area, A 
similar and very widely accepted definition based on the same 
principle, was stated by Holzworth (1967). He contended that the 
mixing layer is that region where surface convection coupled with 
turbulence causes strong vertical mixing. Others consider mixing 
depth or height to be the level at which parcels of air from the 
surface reach equilibrium with their surroundings. Workers in the 
field of large-scale air pollution studies prefer to consider the 
mixed layer and the PBL as one and the same (Laikhtman 1961, Hanna 


1969, Agnew and Jarvis 1972, Munn 1975). 
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From the previous section, the mixing depth has been 
generally defined as the height to which surface air parcels will 
rise. This creates a layer containing most of the surface emissions 
from a wide variety of man-made and natural sources. The extent to 
which the air is mixed becomes critical when air pollution concentra- 
tions are being determined. As the depth of this layer varies, the 
efficiency of diffusion-dispersion processes also varies. Therefore, 
the purpose of forecasting or in some manner determining the height 


of the mixed layer is to aid in the prediction of pollution episodes. 


To2 ) ptructire of ithe Mixing Layer 


The growth of a mixed layer depends on processes which 
are time dependent, and, therefore, the development of this layer 
is also a time-dependent phenomenon. The strength of the mixing 
in the lower layers and the degree of contrast between the upper 
more stable airmass with that below, will determine the intensity of 
the inversion which delineates the top of the mixed layer. In 
theoretical models, Carson (1973) and others, describe the lower 
atmosphere by using three distinct regions (see Fig. 1.1). A 
superadiabatic layer dominated by forced convection develops close 
to the ground. Next is an adiabatic or free-convection layer which 
is capped by a deep non-turbulent zone. Carson's model contains 
a fourth layer which is not as distinct as the others. This is 
the interfacial entrainment layer which separates the free-convection 
region from the more stable layer aloft. This region of entrainment 


varies both in depth and character in the real atmosphere but does 
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act as a boundary to the free-convection layer. 

Carson's model accurately depicts all but the lowest layer 
of the atmosphere as observed by Holmgren et al (1973) who studied 
acoustic sounding and balloon ascent records over Fairbanks, Alaska. 
The proximity of forest vegetation appeared to hamper the development 
of a superadiabatic layer next to the ground. This discrepancy was 
attributed by Holmgren et al, to a warming of the air by the nearby 
forested areas. Carson's model atmosphere describes the real 
atmosphere well, with some modification to the lowest layer needed 
under particular conditions. It will be used as a basis for the 


rest of this study. 


1-3) Review OF Mixing Height Models 


The first attempts to predict the height of the PBL relied 
totally on the theory provided by classical fluid dynamics (Hanna, 
1969). When this theory was applied to atmospheric problems, such as 
calculating the height of the Ekman Layer, estimates were found to be 


proportional to; 


rh 
ae mt (2K/£) * (1.31) 
where Hy is the height at which the wind direction in the PBL first 


coincides with that of the geostrophic wind. The eddy viscosity is 
given by K, mt = 3.1416 radians and f is the coriolis parameter. 
This led Blackadar (1962) and others to propose that Hp could be 
calculated under adiabatic, steady-state, barotropic conditions 


using the formula; 
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Hy = H(US/£))— 0.2 ee (1352) 


which incorporates the friction velocity near the ground, Wee 
with the coriolis parameter f£.. Equation (1.3.2) has a similar 


form to one proposed by Rossby and Montgomery (1935); 
H(U,/£546) = 0.021 (U,/£)sin(a,) (12323) 


where a is the geostrophic wind at the top of the PBL and ay is 
the angle between the direction of the surface wind and that of 
Uy. Equation (1.3.3) was also applied in the following form to stable 


conditions; 


5 AQ] +4 
H(Ug,A0,09) = 0.38 Ussin(ag)/ Aa ? (15954) 


The lapse rate of potential temperature, AO/AZ, and the mean temper- 
ature of the layer, T, are combined with the acceleration of gravity, 
g, to represent the influence of buoyancy forces. Equations Gir 32) 
and (1.3.4) were tested by Hanna (1969). The results exhibited a 
large degree of scatter when compared with measurements of the 
mixing height. In a later study by Carson (1971), the formulae 

were again tested on the same data with similarily poor results. 
Carson attributed the wide scatter to the use of equations based on 
steady-state similarity theory and indicated that time-dependent 
models were absolutely necessary. A similar conclusion was reached 
by Deardorff (1972). The time-dependent equation for the height 


of the mixed layer, as derived by Carson (1973), is given by; 
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h@(t) = h*(t.)exp{28(t, - t)} + 


c 


+ 2exp (~28t) f exp (281) { (H(0,7) - 2u(h,t))/pey(t) fdr (1.3.45) 


c 
Oo 


In- this formulation,._the height of the PBL.-h(t), is related to the 
initial height h(t )) by integration of the change in heat flux at 
Ehe surface,vH(0,7).tosthat,atythe heightah,yH(h,z) for a given) time; 
t. Also included is a constant convergence parameter, 8, describing 
large scale convergence which leads to subsidence. The B- term is 


related to stability by; 
y(t) = (0)exp (BE) (i206) 


where y (t) is the lapse rate of potential temperature and y(0O) is 
the same quantity at t = 0. The net effect of 8 is to increase 
stability due to large-scale subsidence. Deardorff, Willis and 
Lilly (1969) assumed no subsidence (8 = 0) and a constant surface 


heat flux (H,) > whichpreduces, (1.3.4) tos 


2H oe 
) 


h2(t) = (1.3.7), 


tee: 
where y(0) is the lapse rate of potential temperature at t = 0. 
The simplification of (1.3.5) to (1.3.7) produced a model which 
demonstrated success in calculating the height of the convectively- 
unstable boundary layer (Carson, 1973). This formulation (1.3.7) 


should be restricted to days with a dry convectively~unstable 
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boundary layer under almost clear~sky, non-advective situations. 

A similar approach was presented by Summers (1965). He 
considered the advection of a column of air across an urban area. 
A method was proposed which would slowly alter the rural lapse rate 
by adding heat to the air column from below. The amount of heat 
added dictates the height to which the dry-adiabatic layer will 
extend. The level at which the dry-adiabatic lapse rate changed 
to the rural lapse rate defined the height of the mixing layer. Using 
the amount of heat available from an urban area, H, the mixing height, 


h, is given by; 


* . 
De HEU COC) (iv 328) 


where L is length of trajectory of the city, U, is the surface wind 
in the rural area and a is the difference between the rural and 

the dry-adiabatic lapse rates. This method is widely accepted for 
urban heat island situations and has been used by Leahey and Friend 
(1971) and by Hage (1972) with minor modifications. The major 
problem with this model is the assumption that all available heat 
comes from the urban heating rate, H, hence ignoring any other 
sources of energy. 

A similar technique to that of Summers (1965), uses the 
morning temperature sounding (1200 GMT) and replaces the urban heat 
rate by the rate of increase of surface temperature. In this method 
the intersection of a dry adiabat through the surface temperature and 
the morning vertical temperature profile defines the height of the 


mixed layer. Hosler (1961), Holzworth (1964 and 1967) and Miller 
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(1967) used this technique both in calculating mixing heights and 

as a forecast tool (Miller 1967). The use of a dry-adiabatic lapse 

rate requires that certain implicit assumptions be made (see Section 
2.3). <A practical problem in using this technique results from the 

requirement for a morning sounding which may not be available for 


the area of interest, 


1.4 ‘Outline of Method Usediin this Study 


The main parameter appearing in most of the previously 
mentioned models is the surface heat flux. To properly describe the 
energy available to the lower atmosphere, the present model uses a 
surface energy wat nee. In determining this balance, consideration is 
given to the effect of changing cloud cover, varying sun-angle and 
surface albedo. The surface heat flux is assumed to be equal to the 
net radiation, neglecting losses to the ground. Once the surface 
heating rate is determined, a technique similar to that of Summers and 
Holzworth is used. 

Assuming that the morning lapse rate is linear with height 
and that a constant lapse rate for the remainder of the day can be 
predicted, the effect of the surface heat flux can be represented by 
a triangular area on a temperature-height graph (see Section 2.2). 

By using a thermodynamic diagram, area and energy are easily trans- 
formed into an equivalent area, the altitude of the triangular area 
is calculated giving the mixing height. This eliminates the need 
for an observed temperature profile and allows mixing layer depths 


to be calculated for areas where such profiles are not available. 
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Most of the data required for the model are used in the 
radiation component. The site location, cloud amount and surface 
albedo are required together with the minimum and maximum temperatures 
for the day of interest. If predicted or observed heights and temp- 
eratures are available for the 850 and 700 hectopascal surfaces at 
0O and 12 G.M.T. a more detailed form of the model can be used. 

The final piece of information is the date (day and month) which 
gives the correct solar declination. 

The resulting model is time~dependent and produces hourly 
values of the mixing height from sunrise to late afternoon. The 
required input data for this model are available from forecast 


weather offices in the area of interest. 


CHAPTER 2 


ENERGY APPROACH TO MIXING DEPTH 


Za Development of a Mixed Layer 


For the development of a mie layer to proceed, processes 
which destroy stability and therefore produce instability must be 
present. The primary forcing agent is thermal or convective turbulence 
generated by the surface heat flux. An upward surface heat flux is 
created by solar insolation or by cold air advection over a warm surface. 
Once the basic turbulent state is established, it can be enhanced by 
the advection of cold air, radiational cooling aloft or by forced con- 
vection. Wind over rough terrain will add to the turbulence through 
enhanced mixing. In the present model, an upward surface heat flux 
produced by solar heating is considered to be the fundamental energy 
source with mixing enhanced by advection being considered only indir- 
ectly. 

The strength of the convective eeeeed nee depends on the 
time of year, the time of day, cloud amount and the type of surface 
covering. This process is therefore highly variable even over very 


short time intervals. For convective turbulence to be effective, the 
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air next to the ground must become warm enough so that buoyancy 

forces will overcome gravity. This leads to an unstable layer in 

the lower levels of the atmosphere and usually results in an inversion- 
capped mixing layer as described in Section 1.2. Surface heating is 
the primary forcing agent which will develop a mixed layer and must be 
the basic parameter in mixing height models. 

A strong correlation exists between mixing and instability, 
hence a process which alters the ambient stability will be important 
in the development of a mixed layer. Such a process is advection, 
both frontal and non-frontal. Cooling or warming in the lower levels 
associated with frontal passages can very quickly alter the stability 
of the lower eee ee The passing of a cold front can cause very 
stable conditions to break down becoming unstable in minutes. Warm- 
air advection usually has a stabilizing effect but can lead to mixing 
if the warm air reaches the surface and combines with the sensible 
heat flux resulting in a warmer but unstable airmass. The advection 
of cold or warm air over the location of interest will usually alter 
the stability and, therefore, will be of great importance when the 
depth of the mixing must be determined. 

Another process involved in the erosion of the stable air 
above an inversion is caused by a combination of wind shear and 
convective turbulence (Carson and Smith, 1974). This process draws 
heat downward into the lower layers from the warmer, more stable 
air aloft. The effects of this 'entrainment' of air through the 
inversion is discussed by Barnum and Rao (1975). However, the 


exact role of this process in the break down of an inversion is not 
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well understood (Turner, 1973). 

In the high arctic, for at least two of the winter months, 
no incoming solar energy is available. This implies that thermally- 
driven processes resulting from surface heating will not occur. With 
no convective turbulence, mixing of the air must result from either 
advection or radiational cooling aloft, or from a combination of the two. 
In figures 2.1 and 2.2, radiosonde ascents from Resolute Bay in the 
N.W.T. show how the combination of the two processes act to provide 
a mixing depth without surface heating. It will be necessary to 
include both advection and radiational processes in any model which. 
is to be used in the arctic in winter. 

The nee of mixing as described above create a homo- 
geneous layer based at the ground. It has been calculated (R. Portelli, 
1967; C.L. Nortonsand’C.B. |Hoidale, 1976) that this mixing layer may 
extend to several kilometers. Portelli (1976), using Holzworth's 
technique, determined monthly mean maximum mixing heights for many 
areas in Canada. A great seasonal variation can be seen in values for 
the northern locations (Teble 2.1). An equally dramatic range in 
mixing heights due to latitudinal variations is demonstrated by 


Portelli's work (lable> 2.1 and Figures 273-227). 


2.2 Energy and Mixing Depth 


Once the sources of energy have been defined it remains to 
devise a method by which the amount of energy can be translated into 
a workable mathematical form for prediction of mixed-layer depths. 
The first step in this manipulation was provided by Godson (1958), 


who derived a relation between the energy provided to the lower 
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QO) == PK COORD. 


where Q, is the available energy, A is the energy per unit mass and 
P is the mean pressure of the layer of air immediately above the 
earth's surface which is affected by the surface heating, Rewriting 


(2,-2.1), and solvine for A .in terms: of "the available energy, Q 


a>? 
oA 
K=O 28 (222 2) 
P 
when A = joule g 
Qy = joule m 


Pp = hectopascal 


Thesronm, Of Avin(2.2.2)) does. not lend atseligwell tovicalculation. 
Transforming to area in metric units on the tephigram as published by 


the Atmospheric Environment Service (A.E.S.) yields; 


Q 
A = 394,9-—4 (2.2.3) 
P 
when A= me 
Oa = joule m 


p= hectopascal 
Equation (2.2.3) permits one to calculate an area in square meters on 
a standard tephigram given the available energy in joules per square 
meter. However, the area itself is not very useful unless the geometric 
shape which the area occupies can be approximated. 
The growth of a mixed layer (see Fig. 2.8) is frequently 


observed to be accompanied by the formation of a shallow superadiabatic 
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(a) Morning Sounding 

(b) Early Stages of Modification 

(c) Maximum Mixing at Maximum Temperature 

(d) Surface Heat Flux Negative - Mixing Reduces 
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layer which becomes more nearly dry-adiabatic with increasing 
height. Assuming that the lapse rate in the mixing layer can be 
adequately described, usually by a dry-adiabat, an energy-area 
technique can then be employed. The technique amounts to approxi-: 
mating the deformation of the lower atmospheric sounding by equating 
the energy to the area of a triangle, Kagawa (1968). From Neiburger 
(1941) and Holzworth (1967 and 1972), assuming that the sounding 

is affected by surface heating only, the depth of the mixed layer 
can be determined by the intersection of the dry-adiabat through 

the maximum surface temperature with the temperature profile aloft 
from the morning sounding (see Fig. 2.9). 

Alternatively, given the area and surface temperature 
excess above the minimum temperature, and assuming a constant temper- 
ature gradient in the mixed layer one can compute the height of the 
mixed layer without knowledge of the initial temperature profile. 

The only new assumption needed is that the initial temperature 
profile is linear. This technique must be abandoned when surface 
heat flux becomes negative, Businger (1973) describes the collapse 

of the dry-adiabat next to the ground near sunset leaving an elevated 
mixing layer with a stable layer developing beneath. For this reason, 


the present model is terminated when the heat flux becomes negative. 


2.3 Adiabatic Assumption 


Minisonde data from Edmonton and Fort McMurray, Alberta, 
appear to support the assumption that a dry-adiabatic lapse rate will 


develop and persist. In numerical models of the PBL, constant potent- 
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ial temperature layers are developed and maintained for the duration 
of the forecast period (Carson, 1973) in order to describe the lower 
layers of the atmosphere. Most studies, however, have been performed 
in areas that exhibit strong surface heating consequently developing 
dry-adiabatic layers reasonably easily. The present study includes 
periods in late fall and early spring when the surface heat flux is 
weak. Therefore the development of the adiabatic layer is retarded 
and only extends a few hundred meters in many cases. This is the 
critical time of the year for pollution problems when mixing layers 
may extend just above the stack heights. 

For many locations in the arctic, solar insolation dis- 
appears for one or two months in the winter. Unstable layers are 
still able to develop (see Fig. 2.1 and 2.2) even when heat fluxes 
at the surface are negative. These mixing layers must be developed 
by advection-entrainment processes as discussed earlier. Unfortunate- 
ly the present model is not capable of calculating entrainment or low 
level advection effects. These situations will not be considered in 


the present study. 
2,4 .Dit ficultdteswmutnasurtace Heat Klux Medel 


The main difficulty with this type of approach is that 
the surface heat flux is a parameter dependent on many other meteor- 
ological processes. Calculation of the heating rate Henrore depends 
on the model's ability to correctly depict properties such as cloud 
cover, surface temperature and albedo as well as atmospheric stability. 


It is very difficult to obtain an accurate measure of the upward sur- 
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face heat flux when the sun angle is low and the albedo and cloud 


amount are large. Care must be taken to recognize when calculated 


values of the flux are nearing the same order of magnitude as the 


model errors. this peoblem will occur in latesfall) and early spring 


and must be considered when results are being analyzed.- 
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CHAPTER 3 
MODEL DEVELOPMENT 


3.1 Shortwave Component 


This component of the model calculates the contribution of 
the shortwave radiation to the total energy balance. Incoming solar 
energy passes through layers of cloud and then reflects at the earth's 


surface. The available energy is determined by the equation; 


RL = (1 - 94) (Q.+ q).—- Lt + Ly Goi) 


where noes the net radiation, Lt and LY are the upward and downward 
longwave radiation components, respectively. The shortwave radiation 
is ‘divided into the direct: radiation, Q, and the diffuse part, iq. 
As will be seen later in Section 3.3, the albedo of the earth's 
surface, a, is a very important parameter in the calculation of the 
radiation balance. 

In this model, direct and diffuse solar radiation are com- 
bined into the available shortwave radiation, Q,. The available 


surface energy on clear days, Q,”, can be calculated by (Goddard, 


1973); 


73) | 


x 
< 
iy f i) 
we 
Wy) 
£ - 
Bd 
‘ 
= 
& 
i 
s 
= * 
L ‘ 
= } 2* 
j 
4 i 
oe » 
a] 4 ‘ ne: 
' 5 Le ' as 
= My - 
4 
a 5 
Fate “me ba Pe 
‘ Ay £ . FS nia & 
> 
= 
Se wt 
om ) q - 
Wart a , > 
[oe TAS oa * 
‘ 
er, as ¥=5) eit 
ive 


Pree! BGEI515 


= 4 


¢ Pant att i 


daatete’ 


“es 
a Vee 4 ; : ah 
A ary wise 
2 - : 
wt tty z ; Bo 5 : 
AGWLAG bay) . 
ee | 
‘ Boxed! peed Bt 
1 
Ks = Pa aes 
i rid tor eters 
r Aer Wy 
o- s <) 
ba co oe a lag 7 
ne = 
4 ) } ; An Val | 
ce Oh i, pata ae 
. . = at «ER bhea §> @ ee a 
a Views! Bal ae a ae SY eae r. is Labbe ue 
7 x 
ss Ow = * * ete » ; rae ig 
6“ ‘ » r : > * q Mi 
i! = nt 4 iz ' coe | ' Een 
eg s 5 pee Oa 
i : ; A ee ? ; 
> NG ee ait a 2 Ser 5 
rine, Bg nae Sexe aia sa 
a 7 
i oe 5 I 
e= Ln) wee ‘= fi rte 
at 3 4 + J! te i 


as aera ied Lp tena ‘wi 


St eee ROE ed +a 


RA Bt Aa u 


at er otal nt a 


a é . a % 


7 
an 


@au-.QaT cos (Z) (3.042) 


Hence, the evaleie pee is related to the solar constant Qo? the 
atmospheric transmission coefficient T, the optical airmass m 
(m= see Z) ‘and the zenteh angle of the sun Z.” Equation (3.1.2) 
represents the entire shortwave model bee Sei to calculate the 
average attenuated intensity of the solar radiation received at the 
ground (Goddard, 1973). 

The formulation shown above contains implicit assumptions 
with respect to the transmission coefficient, T. Equation (3.1.2) 
contains the simplest and crudest measure of atmospheric transmiss- 
ivity (Stagg, 1950). No attempt is made to isolate particular 
wavelengths or to distinguish T from the true monochromatic trans- 
missivity q*. It should be noted that q~ is independent of m 
whereas T is not, as observed by Stagg (1950) and Latimer (1974). 
These characterisitics demand the calculation of local values of 
transmissivity in order to improve the simulation of atmospheric 
conditions over the area of operation of the model. To this end, 


values of T were obtained by rewriting (3.1.2) in terms of T; 


f=" exp “log pees (Que 3) 


O 


The values for QF were obtained from the Monthly Radiation 
Summary for Stony Plain, Alberta with care being taken to insure that 


clear days only were chosen. After the analysis was performed, a 
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diurnal aasenadade of T on m was observed. However, such dependence 
was considered small enough to ignore when compared with the hourly 
fluctuations. This implies that T is independent of m in this model. 
The main variation of the transmission coefficient was seasonal. Three 
values of T were obtained from the data by averaging over three monthly 
periods. Further adjustment was made by running the radiation model 
and comparing shortwave radiation value (actual and observed). The 
characteristic values and standard deviations are given in Table 3.1. 


Similar trends were observed by Stagg (1950) and Latimer (1974), 


able 3.1 


Seasonal Transmission Coefficients for Northern Alberta 


Season 


Winter 
(Nov-Feb) 


Transmission 
Coefficient 


Spring & Fall 
(Mar-Apr) - (Sept-Oct) 


Summer 
(May-Aug) 


It ayewtla be noted that Q,° calculateduby, (371% 2) e1s 
the sum of the direct plus diffuse radiation for clear days only. 
Modification is required for days with cloud cover. In order to 
simulate multi-layered cloud cover, a model proposed by Pandolfo, 
Cooley and Newburg (1963) was incorporated. The original formulation 


accepted four cloud layers of different height and amount. In the 


29 


- 


tet 
7 ae 
ae Ve 


eh Taras) 


: prey ath? a “* 


= 
Be 
- wg a ih 
Bey e. = i) wi 
: ese : ee PF essa) 
/ : ere ye wey De ore Ane OG 
gig an ae vant hs | 5 ECT ARR rch Gs 
a ae eee 
= ‘ L , = — J a 
Papi ore © sii & ied bao 
ae sSehin re RL oO) Penh 
< (i: 
a ~ ; ee, 5 
14, 4anto tee eet 
4 ys] es - ade 
we a r 2 
vs ' es 
EPP Th — ope 
r - PAL Pilon. 
. ‘t's eee 
~ ie 
4 Se 
¢ >a eer wees hf 
re < ra ~ Hi 
+. > * vies 
% WE si 
1 
is 
& : > 
“eS 
- Te’ ae 
q 
q 3 ol a ore nave » is ae 
=e 1 D+? 
\ + 7 ’ 
’ 
2 eee ¥ 4 
- i; | 
7 | 
‘ 2 ese") ' 
, cay 
yer i 
« a ca ‘ - " , 
4 2 a ss + 
=? 9 4 ; 
ie ye ' i : iy 
me J $¢ f] we 
oer, Soe ae 
= i 
\ C : P | 
1 carla Eis - ee 2 er) ee eee 
‘ s ae | 3 ; 
e * ‘ ha j 
ine: 
rad | 
4 x rs iM 
hao : ah ° ‘ ' a » ih 
ae Pa 3 7 A eee r 
¥ s =) > a en 
lois ‘ 
- Sa a Oe ax = dah ee 2 
Rae yeh geek: 22 ¥ai cae. tad 
Fs : - ei 
; 4 
Pa ey a eae Col sige aes, rea Hy 
$ 
; ‘4 ; 4) a 
2 are Os ae Ce ae ; af aerae 
' é 
¥. ‘ a 


Go 
© 


present moda only three layers are used. The heights of the various 
layers are incorporated by the use of transmission coefficients which 
take on values indicative of the respective cloud layers. Values 

used by Pandolfo et al. (1963), for transmission of shortwave radia- 
ELON sare extrapolated from Fritz (1951) \Gee8lable 1332). These 
transmissivities are in good agreement with those presented by Reynolds 
et al (1975) and Liou (1976) where absorption of shortwave radiation 
by the cloud is also considered. The available shortwave radiation 


reaching the earth's surface, is, therefore, described by; 
Qa = (1 ~ Ny (1-91) 1 = Ng (1-05) 1 ~ N3 (1-3) 0 


or 


Cl - as bE 


| 


Ny is the cloud amount in the ith layer and jy is the transmissivity 


GrEs4) 


for the ith layer. 
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Cloud Transmissivities 


Transmissivities(%) 


Stratus 


Stratocumulus 


Altocumulus 


Cirrus 


The eee in Table 3.2 are taken from Fritz (1951) and from the 
Smithsonian Meteorological Tables, 

The final modification needed to complete the shortwave 
component is the depletion of radiation due to reflection at the earth's 
surface. Values of surface albedo used in the model are those obtained 
from the Smithsonian Meteorological Tables. Further discussion of sur- 
face albedo is deferred to Section 3.3. The available shortwave rad- 
iation in the model is calculated by; 


3 
Qa = (1 - a) a - N;(1-,)) TMcos(Z)Q. SAS 
1= 


3.2 Net Longwave Component 


To complete the surface energy balance, the emission of 
longwave radiation must be studied. The surface of the earth emits 


energy as a grey body following the Stefan-Boltzman law; 


Lt = ect BE2e)) 


where ¢ is the surface emissivity, o is the Stefan-Boltzman constant 
and T is the absolute temperature. The value of the emissivity is 
usually between 1.0 and .85 (Petzold, 1974) with temperature measured 
at the ground level. Considering the operational nature of this 
model, ground temperature was replaced by screen temperature. This 
assumption was also made by Petzold (1974). However, such a crude 
approximation would not be appropriate for microscale modelling. 
Longwave radiation emitted downward from the clouds acts 


as a partial balance to the surface emission. Brunt (1932) and 
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others, modelled downward longwave radiation using both température 

and vapour pressure in their prediction equations. Swinbank (1963) 

was able to show that a more general and equally acceptable formulation 
could be obtained by using a modification of the Stefan-Boltzman law. 
Swinbank (1963) then performed a regression analysis for clear nights 


which gave the following relation; 


(32a) 


where T is temperature in degrees Kelvin. Temperature in (3.2.2) 

should be the temperature of the emitting body. This, however, is 

often difficult to ascertain. Using the theory of longwave trans- 

mission and emission it is possible to show that effectively all of 

the downward flux originates at a level very close to screen height 

(Swinbank, 1963). Hence screen temperature can be used in (3.2.2). 
Combining (3.2.1) and (3,2.2) results in a relationship 

for the net longwave radiation (Li): It must be noted that the down- 

ward flux calculated, as in (3.2.2), is for clear skies only and 

must therefore be modified for cloud cover. The complete longwave 


radiation balance can then be written in the following form; 


Ly = ~eoT? + nk5,31 x 10°14 Te RAS) 


The Swinbank formulation (3.2.2) is, therefore, modified by the cloud 
amount n and by a factor k which depends both on cloud amount and type 
(Reuter, 1951). In a similar formulation by Goddard (1973), the upward 


‘flux rather than the downward flux is modified by cloud amount. Comments 
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by Goddard (1973) indicate good results for Pisa sky situations. 
Cloudy occasions were not tested adequately. 

The method used to calculate the net longwave radiation in 
the present model differs from both of the above methods. An analysis 
of the above techniques indicated that the combination of upward and 
downward fluxes did not correlate as well with L. as did the screen 
temperature. This led to a regression equation for the longwave net 


radiation in the form; 
(B21) 


A similar technique used by Petzold (1974) gave very good results over 
snow-covered surfaces. ‘This approach follows logically after a care- 
ful analysis of the oc ieiilas eo (3.2.3). For any given time period 
less than a day, the only variables are cloud amount and type, and 
screen temperature. Since screen temperature also reflects variability 
in cloud cover, the only truly independent meteorological parameter 

is temperature. 

The use of temperature in (3.2.4) must be investigated 
carefully. For each month of the year, screen temperatures will 
exhibit different relationships with cloud cover due to sun angle as 
well as ground cover. As an example, a temperature of -15C could be 
observed in January under a heavy overcast deck of stratus cloud. The 
same deck of cloud may be accompanied by a 0 C temperature in March, 
with snow-covered ground in both cases. To help reduce this type of 
error the constants in (3.2.4) were derived separately each month. 


The method consisted of determining the constants C, and 
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Cy for the months used in the model verification are presented in 
Table 3.3. The 95 percent confidence limits shown in Table 3.3, 
indicate that the largest variation occurs in March. This is due to 
the difficulty in determining a representative albedo value under 
ablation conditions (see Sections 3.3 and 4.2.4), The error analysis 


presented in Section 4.2.2 confirms that these values are acceptable 


for use in the net longwave component of the model, 
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Monthly Regression Coefficients for 
Net Longwave Radiation Component 


Regression Coefficients 
and 95% Confidence Limits 


March 
August 
September 


November 


The final form of the equation used in the model consists 
of the sum of (3.1.5) and (3.2.4) which equals the net radiation 


(RD? available to heat the lower atmosphere. 


(3625p) 


m 
Ro = (1 = a) (Lik, 5 N,v) 7 cos(Z)Q. + c,T. + Co 
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3.3 Surface Albedo 


Perhaps one of the most important parameters in the radia- 
tion model is the value given to the albedo of the earth's surface. 
The albedo value controls the extent to which the shortwave radiation 
influences the overall energy balance. If this value is in error, 
by more than 10 percent, especially over snow cover, it will drastically 
affect the model output (see Section 4.2). 

The most representative albedo values for ae surface 
types were obtained by a verification procedure using published 
values as a guide. A complete list of albedo values is available 
from the Smithsonian Meteorological Tables. However, most of the 
values are reported as ranges requiring the user to extract the part- 
icular value which best describes the situation being studied. eal 
trials were needed using different albedoes to determine the value 
that gave the best reproduction of the observed net radiation under 
the same meteorological conditions. The albedo values in Table 3.4 
were selected by this method. 

In Table 3.4; several different values are given for snow 
surfaces. Snow has the peculiar property of changing its albedo 
sharply from day to day especially during an ablation period (Arai, 
1966). Values can vary by as much as 50 percent over the space of two 
or three days. A trial-and-error process similar to that described 
above was used to obtain the values for snow cover. Some further 
discussion on their use in particular meteorological situations is 
worthwhile, The value of 87 percent is reserved for days where 


measureable amounts of snow are falling, whereas 81 percent is used 
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Lables 4 


Surface Albedo Values 


Grass 


Albedo (%) 


Frozen Ground 


Wet Grass 


Grass (Trace-Snow) 


Snow (Melting) 


Snow (Old) 


Snow (New) 


Tops of Pine 


Lops oL Fir 
for days where only trace amounts have fallen or days immediately 
following a significant snowfall. When more than two days have 
passed with no new snow and with temperatures remaining below 0°C, 

a value of 70 percent is recommended. A value of 46 percent is used 
to describe melting snow with ‘caper tee well above freezing during 
the peak insolation period. Often these periods are accompanied by 
light rain and snow or rain showers. For,surfaces with trace amounts 
of snow and bare ground visible, 25 percent is the best value for 

the albedo. These were the values used in the longwave radiation 
regression analysis described asceuten 3,2. ‘For othersurfaces), 
the albedo values in the Smithsonian Meteorological Tables would be 


more than sufficient for this model, 
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The preceding discussion of albedo values is primarily 
concerned with nearly flat surfaces of homogeneous cover. Ina 
woodland area, the albedo of the treed and open sections must be 
considered separately. For the present radiation model, the value of 
albedo at the top of the forest canopy is used. This requires the 
calculation of a weighted albedo taking the percentages of treed to 
open areas in the region of interest. A procedure similar to this 
was used in the Fort McMurray study and will be discussed further in 


Section 4.3. 
3.4 Diurnal Temperature-Wave 


An integral part of the radiation model, particularly the 
net longwave radiation vemponanie, 6 an hourly predicted temperature. 
This requires a forecast equation which will correctly describe the 
diurnal temperature-wave. From studies of the diurnal temperature 
trend of Edmonton and Calgary, by Hage and Longley (1967), it is 
apparent that a modified sinusoidal wave would be appropriate. 

Before the diurnal temperature-wave model can be formulated, 
two further assumptions are required in order to specify the times of 
maximum and minimum temperatures. It was found that sunrise correlated 
very well with the occurrence of minimum temperatures except in winter. 
From a study of daily temperature trends for Edmonton, Alberta, the 
months from November through to February require that the time of 
minimum temperature occurrence be predicted. Using a similar study, 
the average time of occurrence of maximum temperature was determined 


to be 1500 local standard time except for the winter months where the 
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average was 1400 local standard time. The only remaining problem was 
to correctly describe the behaviour of the diurnal cycle immediately 
after the maximum temperature was reached. From Businger (1973) 
(see Section 2.1), the mixing layer continues to receive energy until 
the surface heat flux becomes negative. Therefore, the diurnal wave 
is modelled in such a way that the value of maximum temperature is 
maintained until the net radiation becomes negative. Due to other 
factors which will be discussed in Section 4.3, the temperature cycle 
is terminated at this point. 

The general form of the temperature-wave model can be 


written in the following manner; 


T = aM + .5Psin(C) Cra) 


The values of a, Mand P are calculated for each day. M is the 
average value of temperature calculated from the morning minimum and 
maximum temperatures for the day and P is the difference between these 
values. Therefore, M gives the appropriate ordinate value for each 
day and P is an adjustment ronthe amplitude of the wave (see Fig. 3.1). 
The distance in time units from the minimum to maximum temperature 


is varied from day to day by C, which is given by; 
C = 4.7 + (w/(TMAX - TSR))t e452) 


The difference (TMAX-TSR) is the time in hours between maximum and 
minimum temperature which is exactly one half the period of oscillation. 
The value of 4.7 is the phase angle in radians when the sine function 


is at its minimum, 7 = 3.1416 radians and t is the time in hours, 
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Diagram of Diurnal Temperature-Wave Component 
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With all the variables defined, the values of a and b can be obtained 
by solving (3.4.1) for the maximum temperature simultaneously with a 
similar equation for the minimum temperature. This completes the 

diurnal temperature-wave model which enables hourly values of temper- 


ature to be provided internally to the main model. 


3.5 Upper-Air Stability Component 


On many occasions in the spring and fall the surface heating 
will be sufficiently strong to destroy the low level inversion and 
increase the mixing depth dramatically. This process is dependent 
both on the energy input to the lower atmosphere and the overall 
stability of the atmosphere. Therefore a technique is required which 
will monitor the cathe a sins between the stability in the PBL with 
tne stability aloft: Theseritical point will be reached whenvthe 
lapse rate is unstable from the ground to some point above the 
inversion level hence allowing increased mixing. 

This component is incorporated into the model by the use 
of forecast heights and temperatures at both 850 and 700 hectopascals. 
By interpolating values for each hour between 1200 (GMT) and 0000 (GMT), 
a complete set of upper-air data is generated for each hour. As the 
heat flux from the surface increases the Ge eoeh the temperature 
at the top of the mixed layer is compared with the upper-air tempera- 
tures at the two levels previously mentioned. When the lapse rate 
between the top of the inversion and either the 850 or 7/00 hectopascal 
level becomes unstable, the mixing height is increased to the top of 


the entire unstable layer. Instability occurs when the lapse rate 
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exceeds 8.5 C ia (mean value of the wet and dry-adiabatic lapse 
rates). It is reasonable to expect that the hourly progression of 
mixing heights may exhibit a sudden increase due to the presence of 
unstable air aloft. An example illustrating the model prediction in 
an observed situation of this type is shown in Fig. 3.2. 

By adopting this component some of the problems of advection 
are eliminated. In Fig. 3.2, cold-air advection has created an unstable 
lapse rate to nearly 850 hectopascals. The model is able to correctly © 
detect this change by using the upper-air stability component. A 
problem that remains unsolved is the case where advection occurs below 
the 850 hectopascal level but does not actually effect the level 
itself. .Warm-air advection, as shown in Fig. 3.3, is not handled 
well by the model and the addition of the above component is of no 
real help. As can be seen by the actual calculation of mixing heights 
(Appendix D), the depth of the mixed layer resulting from surface 
heating alone will not be sufficient to produce mixing heights of the 


magnitude of those observed, 
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CHAPTER 4 
ERROR ANALYSES AND RESULTS 


4.1 Operation of the Model 


In Chapter 3, development of the components of the model 
was discussed. The coupling of the various segments into a workable 
unit will now be considered. As the net radiation becomes positive 
(after sunrise), the model begins the calculation of mixing heights. 
The energy continues to be.summed and transformed into an equivalent 
triangular area as described in Section 2.2. This transformation is 


explicated mathematically by; 


Rg 


t 
— f(t) dt (4,101) 
ona cS 


A(t) =C 


The sensible heat flux, H(t), is integrated from the time of minimum 
temperature, ty. £0 thestime sol interest, jb. )ihegvalue ior P is 
defined in Section 2.2 and is calculated by assuming a hydrostatic 
atmosphere. The constant C transforms units from specific energy 
to area and has the magnitude 2.53 x 104 ke joule -. 

The calculation of the mixing height, h(t), is a problem 


in basic geometry. A simple formula relating the area to the base 
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and altitude of the triangle would be sufficient but this does not 
allow the required flexibility. Considering that the Beno ere is 
not always dry-adiabatic in ne lower layers, it would be highly 
desirable to have the lapse rate used in the model depend only on 

the particular meteorological conditions present on the day of 
interest. This can be accomplished by specifying the angle on a 
thermodynamic diagram between the base level and the chosen lapse 
rate (see Fig. 4.1). Therefore, the formula used for calculating 

the mixing height, developed by the surface heat flux, is of the fol- 


lowing form; 


4A(t) 
DTsinég 


2A(t) 1 


5 
aera }DTcos6| “sind (4.1.2) 


Bee) = 1 + pT? - { 


As tilustrated in Vig. 4.1, the angle of imterest-is#0. A value 

of 0 = 40 degrees corresponds to the dry-adiabatic lapse rate. The 
value of mixing depth printed out by the model is given by (4.1.2) 
until the entire lower layer becomes unstable (see Section 3.5). 

The program requires the input of basic meteorological 
parameters and geographical information for the location being 
analyzed (see Appendix B). With these values and the model components, 
a value of mixing height is calculated for each hour from sunrise 
until the surface heat flux diminishes to zero. The output includes 
a value for net radiation and solar radiation in joules per square 
meter per hour as well as the hourly mixing height in meters. Pre- 
liminary information in the printout includes the geographical 
location, date, albedo and the maximum and minimum temperatures for 


the location of interest (see Appendix A). The addition of these 


iT i 
‘ ; Ar 
} = ° 
‘ vt Ais) Pe ON RS 
: py : 
a ( ;? 9 ~ t min 
= A ’ : 
ol) COR, Re TUL et ee nekeae 
Ps ; ra Sse aT sash me 
iy * , 4 r % 
i 4 
Vv : aie ¥ a { a! 
t : i + 
‘ . YT VN 
rad 
oY) i ve? A 
; ~ 
3 1 
on i 
~ ari ane 
Ge AUF rR 
A. j ' j : 
ys ) us pig Ot ‘ ~ 1 
© fi b ‘3 y , © 
% ) Be ae pt ty rT .% vu 
. heh aia Si! oi | iia kap 8) kaa 
) a Se : pe 7 | rH 
= Lah, POAC la : 
ee ples 8) I wl F ye lide 32 eek shh ei adel * e 
a we ~ SM + -. ae Seca x he ‘ on ames 1 ie ts se * or wiki AN y maa 
i bf : 3 \ ‘eee ay fs ype erie Ty t ay é 
ht \ 
c = ‘ i i ) 
q oc: eon ee 12 ; ea, ahi i if a “> me + 
S$ : i. a 
= : “i: ; i 4 : 
ur tar git? Vika ee ot shaadi ay Aire i 1 4 
q ath ea Sfi ¥ <a | ie Bing < a oe oe oe. ; 
* j : : - ' 
58 x tae. ee ae | af he) an Ty es as | ay - rey 
| - anny ee | y 
: ie } - ‘ 
: oY ‘ 6 
OT ee ae 
ry ad - ; ' JS aay a a ; 
a = as 7 Al ¥ a 7 : g , ie 
i —_ f Lv ty y re ' (i ; - fi 
‘ f : : arr oy ij i ; ae : 
a i 


Height 


Morning Lapse Rate 


Predicted 


Lapse Rate . 


Temperature 


Pio ae ok Geometry for the Mixing Height Model 
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parameters in the printout will provide a check on the performance 


of the model. 
“4.2. Errors. an the Model 


In the development of the model, approximations of actual 
atmospheric processes created errors. Aside from developmental 
errors there are 'separation errors' introduced by using data from 
different locations within the same region. For example, the radia- 
tion values were taken from Stony Plain, Alberta, whereas cloud and 
temperature data were from Nisku, Alberta (Edmonton International 
Airport). The latter type of error is very difficult to estimate 
but must be considered when assessing the overall inaccuracies of 
the model. 

4.2.1 Errors in the Shortwave Component 

As described in Section 3.1, the shortwave component 
incorporates a forecast of cloud amount, derived turbidity values 
and a modelled solar flux. The errors could be rather large if care 
is not taken in the determination of these components. An error 
analysis was performed using observed cloud cover and temperature 
from Nisku, Alberta with Hadiation data from Stony Plain, Alberta 
used as a verification. 

Assuming that the radiation values would be representative 
of conditions observed at Nisku, this analysis was a test of the 
accuracy of both the cloud model and the solar flux model (Section 
3.1). The radiation model was used to analyze four randomly chosen 


days from each of six months. The months used encompassed each season 
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with an equal number of days in eee A graphical representation 
(see Fig. 4.2) gives the variation of the standard deviation versus 
local standard time (LST). Low error values in the morning and 
afternoon with a slight peak near maximum insolation, 1300 (LST), 
correlates well with the variation of the solar flux reaching the 
ground. The shape of the error curve is therefore considered to be 
the result of larger radiation values giving a higher probability. 
of having greater errors. The mean value of the shortwave radiation 
calculated over this data set was found to be 100 x 104 joules m2 
hr7! with a mean standard deviation of 24 x Wor joules m7 2hr7t (see. 
Rigw 42) - 
4.2.2 Errors.in the Net Longwave Component 

The net longwave component was derived using a least- 
squares regression technique (see Section 3.2). By using the regress- 
ion coefficients with their 95% confidence limits, as given in Table 
3.3, upper and lower bounds were calculated using (3.2.4) for the 
net longwave radiation (see Table 4.1). 

Table 4.1 


Error Analysis of Net Longwave Radiation 


Net Longwave Radiation 
(Joules m72 hr7! x 104) 


March 
August 
September 


November 
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Standard Deviation 


40.0 


(Joules m? x 10° ) 


Rte. 42 


Time (hr) 


Standard Deviation of Shortwave Radiation 
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A value BE 10C was used in (3.2.4) to eliminate the effects of errors 
in the diurnal temperature-wave component of the model. From Table 
4.1, the error analysis indicates relatively good results for all 
months except March. The larger percentage error indicated for the 
March calculation was caused by the rapidly varying albedo which 
entered into the original regression (see Section 3.2). Some of the 
problems that arise in deducing surface albedo during ablation per- 
iods have been discussed previously (Section 3.3). 

4.2.3 Errors in the Net Radiation Component 

In order to determine the error in this component, values 

of the net radiation calculated by the present model were compared 
with observed values recorded at Stony Plain, Alberta. Data used in 
this study excluded all days with snow cover in order to eliminate 
the problems with albedo values. The mean net radiation obtained 


Ze 


using this data set was 34 x 104 joules m “hr ~ with a mean standard 


Aeon (see Fig. 4.3). The standard 


deviation of 12 x 10° joules m- 
deviation of the net radiation did not exhibit a diurnal time- 
dependence. An error of aes magnitude (12 x 104 joules m77hr71) 

was found to cause a 200 meter error in the mixing depth calculated 
in oe ae morning hours. The inaccuracy of the mixing-layer 
height drops sharply with time to approximately 20 meters by mid-day. 
It is interesting to note that the error in the net radiation is 
approximately one half that of the shortwave component. This 
demonstrates the effect of albedo as an ‘error filter' in reducing 


the amount of variation originating in the shortwave component by 


exactly one minus the albedo value. 
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4.2.4 Albedo Errors 

In the preceding discussion, albedo was held nearly constant, 
Unfortunately this constancy of albedo does not occur in actual practice. 
A rather thorough explanation of albedo and its variability was given 
in Section 3.3. It was evident that the required values ranged from 
15 percent to nearly 90 percent. The error analysis, therefore, 
considers three representative albedo values and demonstrates the 
expected error in maximum mixing height that would be incurred by a 
predetermined error in albedo value. Albedoes in error by ten, fifteen 
and twenty percent were used to calculate mixing heights from the same 
data. By using the same data, errors from the other variables were 
eliminated. In Table 4,2, percentage errors are given for albedoes 


of twenty, forty-six and eighty-one percent. 


Table 4.2 
Analysis of the Effects of Errors in Surface Albedo 


Error in Albedo Error in Maximum 
Value (%) Mixing Height (%) 


Albedo Value 
(%) 


As indicated in Section 3.3, the value of the albedo used 


over snow cover (81% range in Table 4.2) is critical in attempting an 
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accurate solution of the mixing-height calculation. This is especially 
significant because an error of 72 percent will make the entire proced- 
ure futile when low mixing heights are considered. By determining 
albedo values, as discussed in Section 3.3, the inaccuracies discussed 
above can be reduced significantly. 
4.2.5 Errors in the Diurnal Temperature-Wave Component 

The diurnal temperature-wave component was verified at the 
same time as the shortwave component. Temperatures used were taken 
from four days in each of the six months chosen from the four different 
seasons. The standard deviation of temperature is plotted versus time 
in Figure 4.4. The peak in the error curve occurring in the mid- 
morning is due to the shape of the model curve (see Section 3.4). 
Although the peak appears rather pronounced, it represents a mesma 
error of only 2.4 C with the average error over the entire period being 
LeAnGs 

The errors in this component could be reduced by the use of 
a more complicated function to model the diurnal temperature-wave. 
However, these errors are considered small enough not to contribute 
significantly to the errors in the net longwave component or the 


model as a whole. 
4.3 Performance of the Model in Northern Alberta 


The performance of the model was assessed by comparing it 
with two simpler models. The simplest approach to forecasting mixing 
heights is to use the temperature profile produced by a numerical 
model such as the Primitive Equation Model (PE). A temperature profile 


from such a numerical model was simulated by using upper air data at 
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1200 (GMT) and 0000 (GMT) which is available from the Monthly Upper- 
Air Summary. Using these data, a five-level temperature profile was 
constructed for those periods required by the study. A second and 
slightly more sophisticated approach, Holzworth's (1964) technique, 
was used after being modified to give hourly mixing heights. The 
surface temperatures required for both models came from the diurnal 
temperature-wave component of the present model. It was hoped that 
this procedure would reduce the effect of errors in the prediction of 
surface temperatures from one of the models to the other. These two 
simpler models were chosen to help ascertain the improvement, if any, 
in forecasting mixing heights using the present model. 

The absolute test of performance is to verify the model 
mixing heights with observed values, Minisonde data were used for 
this part of the study. The main criterion for establishing the 
height of the mixed layer is to determine the level at which the 


stratification changes from unstable to stable. In order to obtain 


a representative lapse rate, the temperature soundings were partitioned 


vertically in equal intervals (approximately 100 meters). This acted 
to smooth the profile. However sufficient detail remained for a 
realistic determination of the mixing depth. The present model was 
used to forecast mixing heights over Edmonton, Alberta for November, 
1974 and over Fort McMurray, Alberta for Neo. Anette and September, 
1976. 

When applying mixing-height models to pollution problems, 
the ability to forecast critical mixing heights is of utmost concern. 


For a mixing height to be critical, it must be within a few hundred 
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meters of the top of the pollution source, effectively trapping the 
effluent in the lower layers of the atmosphere, Under winter con- 
ditions, the mixing heights are often less than the average stack 
height (Portelli, 1975). In the spring and fall, however, mixing 
heights are slightly higher increasing the possibility of critical 
mixing heights. 

The days._used for verification purposes did not include 
periods in which continuous precipitation was falling. Precipitation 
is a problem when a dry-adiabatic lapse rate is being assumed as in 
Holzworth's (1964) technique. Advection days were included in the 
study in order to determine the effects that this fairly common 
phenomenon would have on the model being tested.: 

4.3.1 Fort McMurray Case Study 

Minisonde data for this study were obtained from the Alberta 
Oil Sands Environmental Research Program (AOSERP) and from Syncrude 
Environmental Group (SEG). The AOSERP data (Appendix E) cover a 
period from the 4th to the 17th of March, 1976. From this data set, 
The Great Canadian Oil Sands (GCOS) site was chosen as the most repre- 
sentative location for the area. The GCOS minisonde was located on 
higher terrain, out of the actual valley and was, therefore, deemed 
better than sites in the valley. The SEG minisonde data (Appendix 
E) was recorded at Mildred Lake which is at approximately the same 
elevation as the GCOS site. The entire area in the vicinity of Fort 
McMurray is classified as Boreal forest with approximately forty 
percent tree cover. 


The March analysis covers a period of snowmelt with its 
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inherent difficulties in determining albedo values. Forest cover 
added to this problem requiring that a weighted albedo value be 
determined using the percentages of tree and untreed land in the 
area, The values of mixing height obtained from the three models, 
for this period, can be found in Appendix D. A contingency table 
Ghab lerars Ca). Cb) (c)) was prepared for each of the three models. 
From Table 4.3 (b), Holzworth's model produced 41 percent perfect 
forecasts which was superior to 29 percent by the present model. The 
numerical model produced 46 percent perfect forecasts, however the © 
remainder of the predictions were very poor compared with the other 
two models (see Table 4.3 (c)). A rather strong bias toward under- 
prediction of mixing heights is shown in Table 4.3 (b) by Holzworth's 
model. The present model was free of bias. However, some scatter 


throughout the various classes is evident (see Table 4.3 (a)). 


Table 4.3(a) 


Contingency Table for the Present Model - March 


Observed Mixing Present Model Mixing Heights (m) 
Heights (m) 
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Table 14.3 (Cb) 


Contingency Table for Holzworth's Model - March 


Observed Mixing Holzworth's Model Mixing Heights(m) 


0-100] 100-300] 300-500] 500-800] 800-1200] above 1200| — 


Heights (m) 


0-100 
100-300 
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500-800 
800-1200 


_above 1200 


Contingency Table for Numerical Model - March 


Numerical Model Mixing Heights (m) 


Observed Mixing 
Pen (m) 


0-100 


100-300 


300-500 


500-800 


800-1200 


above 1200 
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Strong surface heating and subsequent convective instability 
dominated the August data period (see Appendix E). As in the previous 
study a weighted albedo was determined but in this case remained con- 
stant throughout the period. Contingency tables (see Table 4.4 (a), 
(b), (c)}) were constructed with the numerical model again being superior 
with 51 percent 'direct hits'. However, the numerical model performed 
very poorly in the remainder of the verification. Holzworth's model 
gave 49 percent correct forecasts which was slightly better than the 
41 percent given by the present model. Bias to under-forecast mixing 
heights was again obvious in Holzworth's model (see Table 4.4 (b)). 
The numerical model demonstrated considerable bias to over-forecast 
heights (see Table 4.4 (c)). No comparable bias was observed in the 


present model (see Table 4.4 (a)). 


Table 4.4(a) 
Contingency Table for Present Model - August 


Observed Mixing Present Model Mixing Height (m) 
Heights (m) 
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Table 4.4(b) 
Contingency Table for Holzworth's Model — August 
Observed Mixing 
Heights (m) 


0-300 


Holzworth's Model Mixing Heights (m) 


above 1200 


300-600 


600-900 


900-1200 


above 1200 


Table 4.4(c) 


Contingency Table for Numerical Model - August 


Numerical Model Mixing Heights (m) 


300-600] 600-900] 900-1200] above 1200] 
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The final set of S.E.G. data used contains the period from 
the 23rd to the 30th of September (see Appendix E). The surface heat- 
ing remained strong enough to 'burn-off' the inversion but the number 
of hours of solar insolation wats beginning to decrease. This created 
lower mixing heights than those observed in August. Fortunately, this 
period was free of measureable amounts of precipitation allowing a 
constant albedo value to be used. As in the other case studies over 
the Fort McMurray area, a weighted albedo value was used. 

Contingency tables have been constructed (see Tables 4.5 
(a), (b), (c)) using the September mixing depths calculated by the 
three models (see Appendix D). Holzworth's technique produced 41 
percent correct forecasts, the numerical model gave 38 percent correct 


forecasts, and the present model gave 33 percent perfect forecasts. 
Table: 4. 5Ca) 
Contingency Table for the Present Model - September 


Observed Mixing} Present Model Mixing Heights (m) 
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Table 4.5(b) 


Contingency Table for Holzworth'd Model - September 


Observed Mixing Holzworth's Model Mixing Heights (m) 
Heights (m) 


0-100 
100-300 
300-500 
500-800 
800-1200 


above 1200 


Table 4. 5c) 


Contingency Table for Numerical Model - September 


Observed Mixing Numerical Model Mixing Heights (m) 


Heights (m) 
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Only the numerical model demonstrated a noticeable bias by consist- 
ently over-forecasting the mixing-layer heights. 
4.3.2 Edmonton Case Study 

Minisonde data were available for a period from November 
19th to December 2nd (see Appendix E) from the A.E.S. Ascents were 
taken at several locations in the vicinity of Edmonton, Alberta both 
in urban and rural sites. Only rural minisonde data were used in 
order to eliminate heat-island effects. The ground was frozen with 
trace amounts of snow. With no appreciable melting or precipitation 
in this period, the albedo value was taken to be constant. As dis-- 
eussed Garlicer, in Section,..4.2,. this isJone of the seasons,.of the 
year when critical mixing heights can be expected and, in fact, this 
was observed (see Appendix D). 

The contingency tables (see Table 4.6 (a), (b), (c)) indicate 
that Holzworth'd model is the best with 75 percent correct forecasts. 
However, the present model and the numerical model achieved 69 percent 
and 63 percent correct forecasts, respectively. It is the bias which 
separates the individual model performances. From Table 4.6(a), no 
bias is shown by the present model. However, Holzworth's model again 
demonstrated a tendency to under-forecast mixing heights (see Table 
4,6(b)). The large values appearing in the lowest class are due to 
the correct prediction of zero mixing heights in the early morning 


hours, 
4.4 Discussion of Verification 


The four case studies covered all periods of the year 
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Table 4.6(a) 
Contingency Table for the Present Model - November 


Present Model Mixing Heights (m) 


| 0-100 | 100-200 | 200-300 


Observed Mixing 
Heights (m) 


above 300. 


0-100 


100-200 


200-300 


above 300 


Table 4,6(b) 


Contingency Table for Holzworth's Model ~- November 


Observed Mixing 
Heights (m) 


Holzworth's Model Mixing Heights(m 


0-100 } 100-200; 200-300 


above 300 
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Table 4.6(c) 
Contingency Table for the Numerical Model - November 


Observed Mixing Numerical Model Mixing Heights(m) 
Heights (m) 


0-100 
100-200 


200-300 


above 300 


except winter. Throughout the verification program very little 
difference was observed between the performance of Holzworth's 
technique and the proposed model. Both models did not perform well 
during the March study. This was caused by problems in determining 
albedo values and the high incidence of advection. The advection of 
cold or warm air effecting the mixing height was observed on 7 out 
of 11 days. The remaining studies were not affected significantly 
by either of the above problems and forecasts were slightly improved, 
The use of contingency tables to display the model perform- 
ances eaten prompt ne deeaaty of accuracy. It must be remembered, 
however, that the classes chosen create arbitrary cut-offs which may 
produce incorrect conclusions in assessments of the oder performances. 
By taking the diagonal elements plus the elements one class from the 


diagonal, an overall judgement can be made which may be more realistic. 
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-Using a contingency table, the significance of the forecasts 
can be checked. For this purpose all of the data from the four studies 
were combined. A contingency table was constructed for the present 
model (see Table 4.7) and for Holzworth's technique (see Table 4.8). 
The significance test is based on the hypothesis that no skill is 
exhibited by the forecasts (the null hypothesis). A new contingency 
table es tan developed using the above hypothesis to derive the 
number of forecasts that would appear in each cell if there was no 
skill. These are the values appearing in brackets in Tables 4.7 
and 4.8. Assuming a chi-squared distribution (Panofsky and Brier, 
1958), the difference between these tables of values can be tested. 
If the two sets Ree are independent the calculated chi- 
square value will be less than the tabulated value for the specified 
significance level and degrees of freedom. The chi-square value for 
Table 4.7 was found to be 133.3 whereas 115.3 was calculated for 
Table 4.8. The value of chi-square with 16 degrees of freedom at 
the 5 percent significance level was found from the appropriate 
table (Mendenhall and Schaeffer, 1973) to be 26.3. Therefore, both 
models exhibit skill. The degree of skill shown cannot be demon- 
strated by this test but can be inferred by considering the various 
diagonal elements. 

The diagonal elements in Table 4.7, running from the top 
left to bottom right, indicate that the present model produced 56 
percent ‘direct hits'. The corresponding value from Table 4.8 
was 62 percent perfect forecasts. If one cell on either side of 
this diagonal is included, values of 89 percent and 81 percent are 
obtained from Table 4.7 and 4.8, respectively. The last piece of 


information that can be gleaned from the tables is the percentage 
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Table 4,7 


Contingency Table for the Present Model - All Data 


Present Model Mixing 
Heights (m) 


0-300 | 300-600 | 600-900 }900-1200 


Observed Mixing 
Heights (m) 


above 1200 


0-300 


300-600 


600-900 


900-1200 


above 1200 
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Table 4.8 


Contingency Table for Holzworth's Model - All Data 


Observed Mixing 
Heights (m) 


See aT, 0-300} 300-600 | 600-900 


0-300 
300-600 
600-900 
900-1200 


above 1200 


Heights (m) 


ih 
(6.9) 


(0.9) 


Holzworth's Model Mixing 
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of unsatisfactory forecasts, more than one class from the diagonal, 
In Tables 4.7 and 4.8, 11 percent and 19 percent, respectively, are 
in this category. Both models gave accurate forecasts between 50 
percent and 60 percent of the time, and gave acceptable forecasts 
more than 80 percent of the time, 

Another aspect of this study is the use of the two previous 
models in pollution episode forecasting. It is the occurrence of the 
critical mixing height discussed earlier which is of greatest concern 
to pollution potential calculation. This critical height will usually 
be from 300 to 600 meters above ground when stacks are considered. 

Low mixing heights where a stack plume can penetrate the inversion 
and situations of unlimited mixing (mixing heights greater than 900 
meters) are also of interest in the prediction of pollution episodes. 

Table 4.9 has been constructed to demonstrate the perform- 
ance of the present model and Holzworth's model in forecasting approp- 


riate mixing heights for the three mixing-height categories. The 


Table 4.9 


Model Performances for Pollution Forecasting 


Mixing Height | Percent Correct Forecasts 


(m) 


Present Model Holzworth's Model 


<300 85 
300-600 40 


>900 58 
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present model performed acceptably in all ranges with Holzworth's 
model demonstrating success only in the first category. This test 
of the models was performed using a very small data sample hence 
firm conclusions regarding the demonstrated skill cannot be made 


with confidence, 
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CHAPTER 5 
SUMMARY AND CONCLUSIONS 


A time-dependent computer model was developed to calculate 
the height of the mixed layer without prior, detailed knowledge of 
the vertical temperature structure. In the development of this model 


three basic assumptions were required: 


(1) the energy necessary to develop the mixed layer 
was derived solely from the surface net radiation 
and was transported upward by convective turbulence, 

(2) the vertical temperature profile from the minimum 
temperature was linear with height. 

(3) the vertical temperature gradient in the mixed layer 


remained constant throughout the forecast period. 


Input for the model included the latitude, longitude and 
elevation of the site under consideration. Surface and cloud albedoes 
were required as were the minimum and maximum temperatures for the 
day. The time of minimum temperature was also necessary for the month 
from November to February. In the remaining months the time of minimum 


temperatures was assumed equal to that of sunrise and was calculated 
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by the model. In order to provide the vertical temperature structure 
aloft, values of the temperature and height of the 850 and 700 hecto- 
pascal levels were included for O500(LST) and 1700(LST). With this 
information as input, the model generated hourly mixing heights for 
the location being studied. 

The importance of the model parameters varied from one time 
of the year to the other. During periods when the ground was snow- 
covered, surface and cloud albedoes became the dominant variables in 
the calculation of mixing depths. In the seasons with no snow cover 
and low albedo values, the maximum temperature for the day was of the 
utmost importance, 

The energy approach modelled the growth of the mixed layer 
ee closely through the summer and atte faviee Ley was round: that 
the detail afforded by the inclusion of upper-air data was needed 
to predict correctly the growth of the mixing layer in the spring 
and early fall. It was, therefore, deemed necessary to include the 
appropriate upper-air data in the model to ensure that the desired 
genes was Sweniiedt 

In order to determine the level of accuracy of the model, 

a comparison was drawn with other techniques, Predictions from 
two other models were verified along with those of the present model. 


The following information was obtained from the test comparisons: 


(1) no significant difference was detected between the 
present model and Holzworth's technique. 
(2) the numerical model gave the worst predictions and 


was strongly biased to over-forecast the mixing-layer 
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heights. 

(3) Holzworth's technique exhibited bias to under- 
forecast mixing heights in all of the case studies 
except September. 

(4) the present model demonstrated no significant bias in 


any season. 


The enhancement of the mixed layer due to cold-air advection 
and/or radiational cooling was found to be a problem an the spring and 
early fall. The inclusion of the 850 and 700 hectopascal heights and 
temperatures partially considered the effects of these processes. 
However, there was no oyun in the results cana 55 cooling 
occurred entirely below the 850 hectopascal level without effecting 
that level. A more detailed vertical temperature profile would be 
required to eliminate this problem. 

The energy approach used in this model has been incorporated 
in many other AOE Te models. However, the present model does 
not require detailed knowledge of an initial temperature profile and 
is, therefore, more broadly Bey tieable than a method such as Holzworth's 
technique. Approaches similar to Carson's model require the measure- 
ment of turbulent perturbations which are often impossible to obtain 
with acceptable accuracy. The present model is more easily applied to 
mixing height determination because it does not require such data. | 
Therefore, the present model is looked upon as a first attempt to 
provide an operational model which is time-dependent and capable of 
producing useful forecasts of the mixing depth using input data that 


are readily available. 
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APPENDIX A 


SAMPLE OF MODEL OUTPUT 


LOCATION 56.DEG. 65.MIN. 111JDEG a 22 IN. 
DATE 447 8 
ALBEDO... 0.19 

= 29. MIN.TEMP.= 15. 


MAX.TEMP. 
TIME NET RADIATION SOLAR RADIATION MIXING HEIGHT 
(LST) (I/O) (Ju? x104 ) (IM) 

4, 84 ~27. De 0. 
5 00 -26. oe 0. 
6.00 -17. 12 0. 
7200 8. 44, 56. 
8.00 14 B52 202. 
9.00 89. 145. 388. 
10.00 120. 184, 549. 
11.00 141. Zaie 690. 
12.00 119. 184. Tide 
13.00 Vig 183. 1098. 
14.00 105. 169. 1098. 
15.00 86. 146, 1098. 
16.00 59. Dei hey athe 
1700 276 138 1152 
18.00 -1. 40. 1151. 
19.00 -24, 11. PU197 
20200 34, 0: 1098. 


MAXIMUM MIXING HEIGHT 1152. (M) CCCURRED AT 17 500(CST) 
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APPENDIX B 


PREPARATION OF DATA FILE 


Description 
File Line or Column Input (all data entries are 
Card Number Format right justified) 
i Laz TZ day of the month 
3-4 12 Month of the yearaClito 12) 
5-7 BO surface albedo value 
8-9 2.0 latitude with degrees and 
10-11 F3.0 minutes input separately 
12-14 oie) longitude with degrees and 
15-16 F2.0 minutes input separately 
17-20 F4.0 leecatvom elevation A. s.La) 
in meters 
21-22 £2 °0 time of minimum temperature 
in hours 
23-25 F3.0 minimum temperature in degrees 
centigrade 
26-28 £3..0 maximum temperature in degrees 
centigrade 
Z 1-3 Lo Sey cloud albedo for first layer- 
first 3) hours 
4~6 Boe2 cloud albedo for second layer- 
first: 8° hours 
7-9 Ei 2 cloud albedo for third layer- 
first 3 hours 
5 1-3 | Roar cloud albedo for first layer- 
second 3 hours 
4-6 Be ee cloud albedo for second layer- 
second 3 hours 
i=9 Powe cloud albedo for third layer- 
second 3 hours 
4 . SAME AS ABOVE for third 3 hours 
5 SAME AS ABOVE for fourth 3 hours 
6 SAME AS ABOVE for fifth 3 hours 
7 1-4 F4.1 850 hectopascal temperature in 
degrees centigrade at 1200(GMT) 
5-8 F4.0 850 hectopascal height (A.S.L.) 


in meters at 1200(GMT) 
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Description 
File Line or Column Input (all data entries are 
Card Number Format right justified) 
8 1-4 F4.1 850 hectopascal temperature in 
degrees centigrade at 0000(GMT) 
5-8 F4.0 850 heetopascal, heieht.CA. Seb.) 
in meters at O000(GMT) 
9 1-4 F4.1 700 hectopascal temperature in 
degrees centigrade at 1200(GMT) 
5-8 F4.0 700 hectopascal height (A.S.L.) 
in meters at 1200(GMT) 
10 1-4 F4.1 700 hectopascal temperature in 
degrees centigrade at O0O00(GMT) 
5-8 F4.0 700 hectopascal height (A.S.L.) 


Example of prepared data file: 


= 


CwomaAnNA MN FWN — 


in meters at 0000(GMT) 


1408018566511122037100015029 


011009003 
009000000 
006018003 
017009000 
009003000 
01621098 
01761098 
00402623 
00522761 
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APPENDIX C 


COMPUTER PROGRAM 


LIST OF VARIABLES USED IN PART 1 


D(I,J)=SOLAR DECLINATION 

BA (I,K) =REGRESSION COEFFICIENTS 
SD(I)=LOCATION PARAMETERS 
AB(I,J3) =CLOUD ALBEDO 

UPT (LT) =GRPER-AIR TEMP. 

UPH(I) =UPPER-AIR HEIGHT 


@! GhiGher ar eve) @@& 


DIMENSION D(12,31),SD (8),RD (12,31) ,TIME(18) ,H (18) ,8(18) , TP (24) , AB ( 
*5,3),PIS (24) ,BA(12,14) ,UPT(4) ,UPH (4) , AH (18) 
COMMON UPT,UPH,SD 
READ (Tate (LJ) pt 112) ,J=1,31) 
1 FORMAT (12 (F5.2)) 
READ. CMO ae OC uem . K) 7; K=1,2) ,12=1,.12) 
2 FORMAT (2F5. 3) 
READ (2,3) NDAY,MONTH,ABO, (SD (L) ,L=1,8) 
3 PORNAT (2U24hoe2 2 F200, F3.0,F2.0,F".0, F220, 273-0) 
WRITE (6,4) (SD(I),I=1,4), NDAY, MONTH, ABO,SD (8) ,SD (7) 

4 FORMAT('1", "LOCATION! ,2X,F3.0,'DEG.',1X,F3.0, *MIN.', 3X,F4.0, "DEG. 
ee pe Oat ENE ey NTs YOY Io, yt TD) 1 ALBEDO N! 2x, Fa C270) 
SeUMAVe TEMP. = 'LP3.0,2%,'NIN. TRAP. = ', P30) 

BEAD (25) ((ab Cusdted= 1,3) ,1=1,5) 

5 FORMAT (3F3.2) 

DOW I=, 4 

READ (2,6) UPT (I), UPH(I) 
6 FORMAT (F4.1,F4.0) 
7 CONTINUE 


LIST OF “VARIABLES USED IN PART 2 


THETA=ANGLE BETWEEN LAPSE RATE AND SUPFACE 
S=SOLAR CONSTANT 

Z=INITIALIZATION HEIGHT 

W=EARTH'S ANGULAR VELOCITY 

F1 AND F2=REGRESSION COEFICIENTS 


C@QIAL@OLe Olle 


THETA=40 .0*0.0174532 

INT=1 

H (1) =1.0 

HF=0. 

NUM=1 

ASUM=0. 

S=487.25 

Z=50.0 

W=.2625 

F1=BA (MONTH, 1) 

F2=BA (MONTH, 2) 

DO 8 IN=1,12 

Po Ss JN=t,31 

RD (IN, JN) =D(IN, JN) *0.0174532 

8 CONTINUE 

G 
C LIST OF VARIABLES USED IN PART 3 
Cc 
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AANQNAAN AQLOQOEOEALOSO Gra 


o) (2) 1@)(@) Leo) (eo) (@) (ei (e) (@) (eel elie) @) (eekie) 


S8=MAXIMUM TEMP IN DEG F. 
S7=MINIMUM TEMP IN DEG F. 
RLAT=LATITUDE OF LOCATION 
ALOG=LONGITUDE OF LOCATION 

HR=HOUR ANGLE OF SUN AT SUNRISE 
TS=LOCAL APPARENT TIME 

TSR=LOCAL STANDARD TIMF AT SUNRISE 
ANP=AMPLITUDE OF DIURNAL TEMP WAVE 
AVG=AVEPAGE VALUF OF TEMP WAVE 


$8=1.8*SD (8) +32. 

S7=1.8*SD (7) +32. 

RLAT=(SD (1) +SD(2) /60.) *0. 0174532 
ALOG=SD(3)+SD(4) /60. 
SM=FLOAT ( (IFIX (ALOG/15.) ) *15) 
TCOR=(ALOG-SM) *4. 

HR=ARCOS (-TAN(RLAT) *TAN (RD (MONTH, NDAY) )) 
TS=12.-57.3*HR*4. /60. 

TSR=TS+TCOR/60. 

LT=I FIX (TSR+.999) 

IP ((MONTH.LE.11) -AND. (MONTH.GE.3)) GO TO 9 
TSR=SD (6) 

LT=I FIX (TSP+ 1.) 

AM P=S8-S7 

AVG= (S8#S7) /2. 

CONST=S8-.5*AMP*SIN (7.9) 


EtST2OF VARLTABLES USED IN PART #@ 


10 


11 


12 


Q=TRUBIDITY 
TMAX=TIME OF MAXIMUM TEMP 


IF ((MONTH.LE.2) .OR. (MONTH.GE.11)) GO TO 10 
IF ((MONTH.LE.5) -OR. (MONTH.GF.9)) GO TO 11 
Q=.76 

TMAX=15. 

GO TO 12 

Q=.92 

TMAX=14. 

GoO'T0" 12 

Q=.875 

TMAX=15. 

CONTINUE 


List OF VARIABLES USED IN PART 5 


APT=LOCAL APPARENT TIMR 
CZ=COSINE OF SOLAR ZENITH ANGL® 

TP (I) =TEMPEPATURE 

SEC=OPTICAL AIPMASS 

PI=SOLAR RADIATION AT TOP OF ATMOS. 
RC=SHORTWAVE RADIATION AT GROUND 

R (NUM) =NET RADIATION AT GROUND 

H (NUM) =MIXING HEIGHT 

TMXH=TIME OF MAXIMUM MIXING HFIGHT 
HMAX=MAXIMUM MIXING HEIGHT 

A=ENERGY-AREA TRANSFORMATION 

ASUM=SUMMATION OF AREA(A) 

DT=TEMP DIFFERENCE (BASE OF ENERGY TRIANGLE) 
AHMX=CORRECTED MIXING HEIGHT FROM UPPFP-AIR TEST 


**MAIN DO-LOOP OF THE MODFL BEGINS KERF** 
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16 
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18 


19 


20 


PA | 


22 
23 


24 


25 


26 


24 


28 


DO 32 I=1,18 
IF (NUM.EQ.1) GO TO 13 
HF=H (NUM-1) 

De(wseT.1) Gosto.(s 

T=TSR 

GO TO 15 

T=FLOAT (LT+I-2) 
APT=T-TCOR/60. 
CZ=(SIN(RLAT) *SIN (RD (MONTH, NDAY) ) -COS (RLAT) *COS (RD (MONTH, NDAY) ) *C 
*OS (W*APT) ) 

LPP. Ghee GO. TO. 16 

TL=1 

GO TO 20 

LE(TOGE. 12.) CO stO 17 
IL=2 

GO TO 20 

TE(T GE. 15-)4 GOnTOnds 
Ti=3 

GO TO 20 

IF(T.GE.18.) GO TO 19 
TL=4 

GO TO 20 

IL=5 

CONTINUE 

TF(P.GE.TMHAX) GO TO 22 
LEMDsCT. A). GOsTG 124 

TP (1) =S7 

GO TO 23 

JST=IFIX (TSR) 
ARG=4.7+3.2/ (TMAX-FLOAT (JST) ) ¥FLOAT (I-1) 
TP (I) =CONST+.5*AMP*SIN (ARG) 
T16=TP (I) 

GO TO 23 

TP (I) =S8 

CONTINUE 

SFC2127C7 

IF (SEC.LE.5.76) GO TO 24 
SEC=5.76 

RI=S*CZ* (Q®*SEC) 

TF (RI.GT.0.),.GO TO 25 
RI=0. 

CONTINUE 
RC=(1.-AB(IL,1)) *(1.-AB(IL, 2) ) * (1--AB (IL, 3)) *RI 
RIS(NUM) =RC 
TEMP=.56* (TP (I) -32.) 

R (NUM) = (1.-ABO) ¥PC+ (F1¥TEMP+F2) *4. 186 
IF (R (NUM) .GT.0.) GO TO 27 
IF (T.LT.TMAX) GO TO 27 

GO TO(26,27) ,INT 
TMXH=T-1. 

HMAX=H (NUM-1) 

INT=INT+1 

CONTINUE 

ZM=(Z+2. *SD(5)) /2. 
PM=1013.25-. 11*2ZM 
A=7.09*R (NOM) /PM 
EP(AsGT.0..). GO. TO 28 
Te(veLT.12.) GO: TO 29 
ASUM=ASUM+A 

PTP=TP (I) 

FSD=S7 

DT=ABS (FTP-FSD) 

IF (DT.FO.0.) GO TO 29 
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29 
30 


31 


32 


SC=(2.*ASUM) / (DT*SIN (THETA) ) 

SA=SOQRT (SC#*24 DT**2-2. ¥SC#DT*COS (THETA) ) 
PHI=ARCOS (( (SA*®*2) + (DT **2) ~ (SC#*2)) /(2.0*SA*DT) ) 
H (NUM) =2104. *SA*SIN (PHI) 

Z=H (NUM) 

IF (H (NUH) .GE.0.) GO TO 30 

H (NUM) =0. 

Z=H (NUM) 

GO TO 30 

H (NUM) =HP 

TIME (NUM)=T 

HMX=H (NUM) 

CALL TEST (T, TEMP, HMX,AHMX) 

IF (AHMX.LE.0.0) GO TO 31 

H (NUM) =AHMX 

CONTINUE 

IF(T.GE.21.) GO TO 33 

NUM=NUM+1 

CONTINUE 


**END OF THE MAIN DO-LOOP** 


33 
34 


ze 
36 


37 


WRITE (6,34) 
FORMAT ('0", 1X,'TIME ',4X,'NET RADIATION ',4X,"SOLAPF RADIATION',4X, 


*OMTKING HEIGAT!' /* *,1X,° (LST) ', 5X," (37M X10 )* 9X," (S78 XOUOM Eat 2X 
Sat) 


JN=NUM-1 

NUH=JN 

DO 36 IT=1,NUM 

WRITE (6,35) TIME (IT), R(IT),RIS(IT) ,H (IT) 

PORMAT (104, 1X,F5-2,6X ,P6-0, 14X,F5.0,14X,F5.0) 

CONTINUE 

WRITE (6,37) HMAX, TMXH 

FORMAT('-"," MAXIMUM MIXING HEIGHT ',¥5.0,'(M)","'OCCURRED AT ', 


*F5.2,' (LST) °) 


STOP 

END 

SUBROUTINE TEST(T,TEMP,H,AH) 
DIMFNSION UPT (4) ,UPH(4) ,SD(3) 
COMMON UPT,UPH,SD 
DLP=-0.0085 

AH=0.0 

RTS= (UPT (2) -UPT(1))/12.0 
RH8= (UPH (2) -UPH (1)) /12.0 
RT7= (UPT (4) -UPT (3)) /12.0 

RH7= (UPH (4) -UPH(3))/12.0 
TOP=TEMP-0.0098*H 
T85=UPT (1) + (RT8* (T-5.0)) 

H85= (UPH (1) + (RH8* (T-5.0)) ) -SD(5) 
IF (H.GE.H85) GO TO 1 

UPLPS= (T85-TOP) / (H85-H) 

IF (UPLPS.GE.DLP) GO TO 2 
AH=H85 

T7 0O=UPT (3) + (RT7* (T-5.0)) 

H7 O= (UPH (3) + (RH7* (T-5. 0) ) ) -SD(5) 
IF (H.GE.H70) GO TO 2 

UPLPS= (T70-T85) / (H70-H85) 

IP (UPLPS.GE.DLP) GO TO.2 
AH=H70 

CONTINUE 

RETURN 

END 
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APPENDIX D 


OBSERVED AND CALCULATED MIXING HEIGHTS 


Values of observed and calculated mixing heights are given 
in Tables D1 - D4. Included in the tables are the results from the 
present model (model A), Holzworth's model (model B) and the numerical 


model (model C) at times corresponding to the observed mixing heights. 


Table D-1l. March Mixing Heights 
Se a ee ee ee ee i) ee ee 


Month: March Location: —PortsMcMurray 


Date Timecesl) 9 Observed Model A Model B Model C 
Heights (m) Heights (m) Heights (m) Heights (m) 
ne a PRE RET oS)! Re A NM OS a 


04 0735 0 0 0 0 
1100 50 260 100 0 
1300 740 311 . 280 IPRS Y 2 
1500 770 1150 310 OT, 
05 0830 0 68 0 0 
0930. 0 Loe 50 1746 
1100 1280 299 310 1754 
~1300 1600 LEL3 490 1765 
1500 1600 1126 900 2500 
06 0930 1200 500 700 1832 
1100 1800 652 1200 2500 
1300 1800 743 1350 2500 
09 0930 120 114 200 0 
1030 200 177 260 0 
1100 0 204 290 2500 
1345 400 310 590 2500 
10 0930 100 200 0 0 
1100 200 310 160 1017 
1300 600 401 330 2500 
1a 0830 0 68 0 0 
0930 200 170 20 0 
1100 500 276 80 547 
1300 1600 378 530 2500 
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Table D-1. 


March Mixing Heights cont'd, 


—_—— eee 


Date 


Ms 


14 


1S 


16 


sgl 


Time (LST) 


0830 
0930 
1058 
1300 
1500 


0930 
1100 
1300 
1500 


0930 
1100 
1300 


0930 
1100 
1300 


0930 
1100 
RS HRS) 


Observed 
Heights (m) 


110 
400 
1100 
1240 
1600 


220 
200 
220 
1600 


90 
140 
1600 


Model A 
Heights (m) 


0 

20 
qZ 
380 
1108 


105 
174 
248 
Ze7 


200 
309 
382 


aie) 
364 
463 


Model B 
Heights (m) 


a 


Model C 
Heights (m) 


0 
2500 
2500 
2500 
2500 


0 
2500 
2500 
2500 
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Table D-2. August Mixing Heights 


Month: August Location: Fort McMurray 


Date Time (LST) Observed _ Model A Model B Model C 
Heights (m) Heights (m) Heights (m) Heights (m) 


ee 


14 0605 0 0 50 0 
0936 297 454 250 >2600 
1304 807 1098 1060 > 2600 
1435 935 1098 1250 >2600 
1802 >2550 1150 1330 > 2600 
15 0605 0 0 30 0 
0938 Die 420 190 0 
1301 680 834 670 >2600 
1800 1700 1201 1020 >2600 
16 0609 0 0 0 0 
0927 977 674 220 >2700 
1254 1105 1193 470 >2700 
1427 7 1600 1375 800 >2700 
1801 >2550 1545 960 >2700 
17 0610 595 0 0 >2600 
0928 807 345 120 >2600 
1255 1105 997 700 > 2600 
1425 >2550 1011 1040 > 2600 
1806 1063 1127 1080 > 2600 
18 0601 - 680 0 0 >2600 
0928 638 367 150 >2600 
1258 680 879 730 > 2600 
1425 1090 1061 790 >2600 
1759 2210 1172 810 >2600 
19 0608 i. 0 0 0 0 
0931 297 673 150 >2600 
1256 2040 1074 520 >2600 
1420 1190 1242 1350 >2600 
1811 1360 1424 1520 > 2600 
20 0613 0 0 0 0 
0943 127, 291 120 >2600: 
dere >2550 1030 650 >2600 
1444 1402 1036 1150 >2600 
1805 >2550 1046 1520 >2600 
21 0616 iy 0 0 0 
0943 297 502 170 >2700 
1316 765 1118 690 2700. 
1441 >2550 1125 1030 >2700 
1801 >2550 1136 1550 >2700 


i 


Siren nnn 


Table D-3. September Mixing Heights 
Month: September Location: Fort McMurray 
Date Time (LSP) Observed Model A Model B Model C 
Heights (m) Heights (m) Heights (m) Heights (m) 
23 0730 0 90 30 0 
1039 213 346 310 1129 
1342 297 450 870 1135 
1427 297 476 1010 1136 
1807 iy | 467 1060 >2700 
24 0725 0 43 60 578 
1043 170 450 350 >2700 
1338 340 558 710 > 2700 
LaSH 425 595 800 >2700 
Lou 850 5795 840 > 2700 
25 0728 0 190 0 0 
1050 255 718 200 576 
1350 22 800 500 578 
1431 275 900 520 > 2700 
1810 977 1024 520 >2700 
26 0716 127 15 0 0 
1056 127 378 340 569 
1355 L275 1130 >1200 >2700 
1440 1063 1130 >1200 >2700 
1807 1387 Ua. >1200 >2700 
ZF 0748 0 75 0 0 
1057 510 318 360 0 
1330 22 738 870 >2700 
1430 850 1081 >1200 >2700 
L7oZ >2550 1081 >1200 >2700 
28 0734 0 45 0 0 
1104 170 B75 500 >2700 
1341 255 875 >1200 >2700 
B35 >2550 1080 >1200 >2700 
1812 >2550 1080 >1200 >2700 
29 0754 765 0 0 >2600 
1109 807 1024 650 > 2600 
1353 1530 1025 >1500 > 2600 
1807 1105 1025 >1500 > 2600 
30 0756 0 153 1180 0 
1102 510 omy 1820 > 2600 
1354 1063 1054 >2000 > 2600 
1445 1400 1055 > 2000 > 2600 
1808 >2550 1060 >2000 >2600 
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Table D-4. November Mixing Heights 


Month: November Location: Edmonton 
Date Time (LST) Observed Model A Model B Model C 


Heights (m) Heights (m) Heights (m) Heights(m) 
ere ee ee ek ty erm en eih Ol) Ee eRe a cee tS Ceo lee eae be Ot 


22 0700 0) 0 0 0) 
0800 0 0 30 0 
0930 36 ay 40 0 
1100 84 174 68 On 
23 0800 0 0 0 0 
0930 0 64 40 0 
1100 44 196 100 0 
1230 116 260 256 0 
1400 260 Zoe 800 Undefined 
26 0800 0 0 0 a0) 
0930 0 86 0 0) 
1100 ; ZZ PAS | 64 270 
1400 436 437 520 275 
Pad 0800 0 0 0 0 
0930 20 18 0 0 
1100 160 68 60 0 
1230 52 127 144 0 
1400 264 147 220 Undefined 
28 0800 0 0 0 0 
0930: 28 45 0 0 
1100 80 147 32 0 
1230 verLaO 187 108 0 
1400 220 205 176 0 
29 0800 0 0 0 0 
0930 90 54 54 0 
1100 124 190 132 0 
1230 150 218 164 0 
1400 168 Parad 228 0 
30 0810 0 0 0 0 
0930 88 68 0 0 
1100 152 247 20 0 
1230 120 232 108 0 
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APPENDIX E 


MINISONDE DATA 


Table E-l. March Data, - Fort: McMurray 
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Table E-1. March Data - Fort McMurray cont'd. 
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Time(LST) 0735 1100 1300 1500 
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Level Height Temp Height Temp Height Temp Height Temp 
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43 1443 Sede, 
44 bay 2 -19.8 
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Table E-2, March Data - Fort McMurray cont'd. 
Time 0830 0930 1100 1300 1500 
(LST) 
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Level Het Temp Het Temp Het Temp Het Temp Het Temp 
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Table-E-3. March Datai- Fort. McMurray 


DATE: March 06/76 Locations..GCOS 
Seer aes cui Ramee dee aCe eae a te eePeeRe nN ee ee ee ee tee ee 
Time(LST) 0930 | 1005) 1300 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
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39 55 -20.0 2107 -17.0 
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Table E-3. March Data - Fort McMurray cont'd. 


Time(LST) 0930 1100 1300 
a 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
a a a a No Sa dl Ms RB di 
Al 1595 =19°0 2197 Pees 
42 1625 =13.7 2243 RSS) 
43 1665 -18.6 2288 pale fey 
44 1685 -18.9 2333 =i oO 
45 ye? =193 2319 SSS 
46 1750 -19.4 2424 eles) 
47 1782 = Ald, 2469 eG 
48 1814 = 0 2515 -18.0 
49 1847 = 20 03 2560 ~18.1 
50 13/9 =20.3 2649 —13.3 
ow ES Tet -20.6 2/38 -18.6 
Sys 1943 AUR) 2827 -18.8 
53 1976 -20.8 2916 erie le 
54 2008 2k 3006 -19.4 
Sy) 2040 -21.0 3095 =19'..6 
56 2073 a Se) 3184 ue Rs} 


Table E-4, March Data - Fort McMurray 
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Se Oe eee te ee ee Se 
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3 433. -19.4 446 ~-16.5 438  -14.6 418 - 9.5 
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29 AB OM eel ae) A341. | LG Ie 29 13.1 9 1329 = 14.9 
30 USS7eH13280 13/1 § Lae e759) 13.58) S75)  b- 1620 
SM PAZ 43 397) | — ESO eh Zoo as Oe eG 2h) OE 
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Table E-4, March Data - Fort McMurray cont'd. 
i ee 


Time (LST) 0930 1030 1100 1345 
EE a a ae Se ee ee 
Level Height Temp Height Temp Height Temp Height Temp 


(M) (C) (M) (C) (M) (C) (M) (C) 


4) 176d 19. 3 1588) -18.8 
42 179549) -19.8 HO2087 —19.4 
43 1829>"=20 74 165i. 19.8 
44 1863) Me2 1.1. LOGSy n= 20.7 
45 Tee 20 18 
46 1746 -21.4 
47 Lee 9 
48 1808 -22.4 
49 US4O009r—=23 7.1 
50 LOvae wn 23.357 
Dt UVOS Aw eH 2G J 1: 
SYS 1934 24.7 
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Table E-5. March Data, - Fort McMurray 
aor ST SE Ne Eo ee ee 


DATE: — March 10/76 Location: GCOS 
Nee ee ee ne 
Time (LST) 0934 1100 1300 
a a re a eee Ne 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
0 368 cdl a) 322 -14.9 320 - 9.4 
Ae 380 -18.2 353 -16.2 349 -10.8 
2 393 -18.8 384 -15.6 a7 9 -11.4 
3} 425 -18.8 415 =o 410 -11.6 
4 456 LER) 447 = OL Gad =i. 9 
5 487 =17.1 474 = 1D 473 -12.0 
6 DLY -16.4 SORE -15.8 505 -12.4 
i, 549 os) 528 -15.6 538 -12.6 
8 580 -15.4 554 =D) 40 Sw) 12.7 
8) 610 -15.6 587 -14.3 602 =13%2 
10 640 -14.6 619 -14.5 633 =13.) 
11 670 -13.3 651 =o 665 -14.0 
q2 700 =i 362 683 -14.2 697 -14.5 
3 731. lb T: 716 -14.6 724 -14.7 
14 763 -14.7 748 cool no yee ou -14.9 
1h) 195 -15.4 778 15a 7 778 alte Gat 
16 826 =15.6 808 16.1 806 <1 2 
17 838 -16.5 833 = 15.17 
18 869 -16.6 860 -16.0 
9 899 -16.6 887 -16.0 
20 930 =16...6 SE Gy 
24 960 -16.6 943 -15.9 
22 990 —-16.7 9735 =—15.3 
2 1021 -16.7 1004 -15.1 
24 1051 -16.4 1036 -15.0 
25 1082 ilies eee! 1068 =15..2 
26 12 -14.4 1099 = yD 
pS) 1142 -14.5 E31 -15.6 
28 : ALS) -14.5 1163 =15.9 
29 1203 -14.4 1194 -15.8 
30 1234 -14.3 1226 -16.2 
Si 1264 -14.5 1258 -15.8 
32 1294 -14.3 1289 -16.1 
33 2S -15.0 to. -16.3 
34 L355 -15.1 2353 -16.7 
35 1386 -15.1 1384 -17.0 
36 | 1416 -15.5 1416 -17.4 
37 1446 -15.9 1448 -17.7 
38 1477 -16.2 1480 -17.6 
39 1507 -16.5 Nog Af -18.3 


40 1538 =1G.9 1543 -18.3 
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Table E-5, March Data - Fort McMurray cont'd, 
ee ee en em See SAN ems Le Ne 


Time(LST) 0934 1100 1300 
PI 1 As a GS RN Fa SS al th a 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
4l 1575 -18.6 
42 1606 -18.8 
43 1638 -18.9 
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Table E-6. March Data — Fort McMurray 
er es eS ee a ee 


DATE: March 11/76 Location: “GCOS 

re Ne ee ee ee ee ee a i OR el ey Se le 
Time(LST) 0826 0931 1100 1300 
A Ea ee a ie sR rd a tal cll tata 


Level Height Temp Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) (M) (C) 


0 B22 aloes S19) =I 6 EVANS Aces) 3200 15). 

a 360 -14.0 362: =—12.2 360" = 9.6 3669) =—"95K.0 

LE SEE Pan is, 404 -12.9 oo ee LO. 2 405 - 5.4 

3 8 1 a mae by ree’ err 30 438° +10.5 LE Grmmtan 5 29 

4 G3, > =43..0 457 -13.3 Ld Mee. 0 491 - 6.4 

3 LON = L205 WoT Cees ca 489 -11.4 9345) = 16.8 

6 fire) Me Ge gk DL6, ~=13.6 5027) =11.8 Gy | 6 

7 Sol" (=13°50 545 -13.8 540 -12.1 6225 =a 

8 Ot ea 13 22 DD.) laa, DOM me 29 6655 ="3.38 

9 O40") wel. 7 604 -13.7 GLO Meta 3.3 1O2Se OZ 
10 682) -13.1 634 -13.6 654 -14.1 150m) 20956 
ae INS 1268 6637) S=13.7 689) (1454 774 ~10.0 
We 149. -=12.8 692. §=13,7 725 lee 810 -10.2 
1 pice, Meese ah 7225 al SiS [De 844 -10.3 
14 816 -13.0 POL pate 2 186° -16,1 Sia eH LO9 
uS 852) -12.5 781 -14.4 816 -16.5 7OSsee—LOn9 
16 884), 2.6 810 =14.7 847 -16.4 OO r= Live. 
iy 92308 -=12.5 840 -14.6 Os / aah Eee 9 Ome Vie 
18 999). 1258 S60.) 5-14.40 DOC lat) LOOM 16 
a9 995 12.7 898 ~-14.6 25 Oe a.0) 9 1029") = 57 
20 LOSOR jo ai361 928 ~-14.8 GOle Lote LODO elias 
yaa 1066: -135.2 957, Lora SIGE 16.4 B LO797 hii. 9 
Ze IO 2 = 13.01 986 /-1554 4) 10277) 16.6) ~11000  .—1270 
23 Ue ee OL VOLS: 1a LOS ee lO Oke 2 2 
24 LIS a1 364 L044. “15.4 1086) 5-17.20) 1142  =12.6 
25 Ome oe LOY 39) Fale ee eer OM Ratt One (22.6 
26 E2405 =14.0 L102) 14.89 1145) '-16.0 (1199> -=2278 
24 VZSOMewam bea: TLS 14.7 Sell / 2 aoe 1226. 12.8 
28 PS LOme-14so $1160) 14,4 2007 me oe 2  ID56  — 12.6 
29 ToS eae) 1189 mA 227 I 16 AE 2B O27, 
30 3G 7h 15.5. L218) ’eE4o2ey elo oe5 11676) el 313 5 9 
34 1423 . -15,9 L247) (je LOZ R289 Wie kOe VSG 2e taped 
32 PGS Oe wea r6.3 . §L276) Slat oe IS 1LOr =i 2 eres eres oe 
33 TA9R TW S=1656>° 1305 | =1428y 1338") 1755) 1399) 2-138 
34 LOSOm@ el: 354 . \-LSe On Is657) Vv Uihoo a a2 oa. 
35 ooo t7i,5- ° 1363. ve Loree eEooS (1779 1461 <-14.5 
36 LGOLe =17..9 1392). =155608) La2d 2 =1S eta o se aS 
oy) NGO mera s.o:, S42) a lao9 L448) i=1875" (elo 2o ai. 2 
38 LGiomee-tOoy), L450, | —1OC4 a tia 60) 919.0 eee O62 eer — 16 
So: POSmera=t950'" L479 ~=L65 6, sal 50 30-) 19 Ore oI 6ee—10., 1 
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Table E-6. March Data - Fort McMurray contild. 


Time(LST) 0826 0931 1100 1300 
Level Height Temp Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) (M) (C) 
41 de 7 -17.2 1663 hdegel. 
42 1566 -17.7 1697 -17.6 
43 7st -18.1 
44 1764 -138.8 
45 1798 -18.2 
46 1832 -19.7 
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Table E-7. March Data - Fort McMurray 
i ee) rem a, 


DATE: March 12/76 Leecation: GCOS 
Time 0830 0930 1058 1300 1500 
(LST) 


SS 


Level Hgt Temp Het Temp Het Temp Het Temp Het Temp 
(M) (Cc) (M4) (cy (M) (Cy) 5 @D (Cc), (4) (C) 


ESSE 


0 S24) iO Oy, O20 =) 625" 319) = 2 Omgeze =42.6 321 = 145 
af S16 Week 505 = 6.9 364 =) Bab eeso “oed.6 ~ 359 = 1.9 
Z EIN Nd SS 7 EM EP me A 0 PR eae WO, SSS San A On 1 Fc ROW 
3 HOS eee G0 S— 7.5. 412) += SiOAOr = 9354 (434 = 2.6 
4 SOD) tee a De at 72 534 “= DON Pa2] y—9 388 PAP (= 2.9 
5 DH eae DOS =i. 3 8596 =) 6a eeeas2 eer4 509. =..3..3 
6 209) Gao eB te 6° 6659" = Pa ee 5 548" =)368 
7 O29 Wai e504 7.6, 108" -) 7.80 B5lG 94S. 582) 6-452 
8 609 (eae 59 St 1. 157 \=<G 4aee > SUtmeen sy 0) OL = aba 
9 109 Ve On 6200 597 59.80) = Be Sao n tw oe Gn.) 652 = 95, 1 
10 JET eWeno OSs 2 8.1 856° =19.3 633° 9=76.3 79687 = 524 
nL TIS een 6S) =) 8. 2 903) = 1940) PO OU O 2 tbs 
a2 O42 = 26. LB 8.2 9547 FLO BO a Orel L TL DE eater Sa 
13 SOCe emorcm ao = Jeo 984 102) asus eo: 3 792) 06 
14 94 79 = 157 1018 2Sk0c0 WiG2, Ba 79 827" — "659 
| Bes 98071 = 75 SHOR io | L057 10 (OPM Cron Oe? I SO2, a2 
16 1026) ee Sar 2864 L085" (a0 A a hop Oe FG Oa O97 a. =a 
L7 NO72, eee alk; = "836 ~ 1107 eS 1LO.3 ¥e9se aoe 6930 fH 769 
18 LEO cmon 902) = 36.9") 1128) 1002) 93d) Ban 9.6 1 962") = 6e3 
19 ALTOS =e eG 929 4952" TI50? =1O 82 0968 \=89.8 92994 5 6.7 
20 L212 O56 2974! 117 S10 83 1005 (9-20, 2) 51026" 2920 
Zi L251 =F SiO 2-968 = "925° 1204) a1l0e5 1042) 1055" 1058) = 954 
22 130342932 1237 =10'6™ 1080 >¥-1059" 1090, =" 9.8 
23 1349" =1001 1270: Sees Se 7 ee eS e225 Oe 
24 1395) =2£0%5 1303" Vela ye hl 5G? eet AS RH LOeS 
Zo 1441 -10.8 ES39 SLi SEL ee Oy eho) Wa LOL.S 
26 1483 -11.1 L370? =1 Gea 28, a2 2d 2 
2] 1534 ~1i4 T40S) Pol S e260" = Po eO) ni 241s atl 
28 1580 ‘=-11.5 TAS 6% R= V2 OS ESI OF e035 20 PIZ695 Te 11s 
29 1626) =li7 1469 SUZ te nS la tS 1293 set? 
30 1672 -12.0 L502 297 81392 150381327 12.4 
a 1718 =1203 1535, Vale eae 2 = Lind, IO, = 1256 
ie 1764 -12.5 ~ 1568 -llv7 1473 -13.4° 1385 ~-13.2 
a8 1G1T 1256 L60) S=147 6 E514 ads 1 42 3) eds 
34 1856 =12.9 1634° Sa G1 555) ee L355) (2460) p13 8 
35 1903. =i5.5 1667 ° Well. 8 84596 @=13.5 41498 1453 
36 1949 1391 1/06, Sel 3 21637 W215. alo relay, 
37 1995" =1352 L734 Fal 2. 7 1G ie Hl SNe loss 
38 2041 - =135.3), L767 iS 620 L716 e139 6m. 4 
oe, 2088 =1355 1800 © =1307" VW 759°=14..00 164 Se-=15 05 
40 2134 -13.6 18337 5=14 J2 


1800 =-14.1 1686 -15,3 
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Table E-7. March Data - Fort McMurray cont'd. 
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Table E-8. March Data - Fort McMurray 
a 


DATE: March 14/76 Location:~ GCOS 


a) 
Time (LST) 0930 1100 1300 1500 


Level Height Temp Height Temp Height Temp Height Temp 


(M) (C) (M) (C) (M) (C) (M) (C) 
ee Ls 


0 324 -16.0 O2he =h23 B24. 4-4 8.6 SOE Boe leusce: 
i 365) wea Low Bod 13.0 SO Aer = 195.5 427 - 7.0 
2 402 --16.4 Sy iS mis 5/72 414° -10.3 545°° = 7.4 
3 44]  -16.4 401 -=13.7 454 -10.8 i ei Ow 
4 480 --16.5 23> '-1358 LD etl LE 608° — 8.6 
5 5065 =16.6 449, -14.1 534° #115 638— = 6.9 
6 535 -16.5 meio .=L433 Somat le 2 6677 oan 
7 568. -16.1 B08.) (L438 606 -11.0 OO aio 10 
8 60M a= 15 53 542 -13.9 639 FO. 7; 42605 = 1908 
9g 633 07ash 4.2 308 —_ -12.6 6/25 1 OKs 756 =-10.1 

10 66655 ~=12.5 pd . -11.6 10> Pelin2 LOD 4 = LO 

11 696" 21.6 6297) §=13.50 133.9) eT 133 OLD g=L0a3 

12 IVES bel de 6521) 1-15 .4 TTA oats 844 -10.6 

13 FOE ses 1D 685.) =13:6 1990 ia, C7358 = LOes/ 

14 TID ee 3 71S) = 14.0 826 ~-12.0 S02 -=—1079 

15 824) =11.5 TDD) (=1453 866 -12.2 EOE ceili ks a 

16 O5ou matte 791 -14.6 JOO a2 99 a LATA 

17 S86" =11.9 813 -14.8 94 2.6 ket) SAS: 

18 SU ig iy O° Saal lear 325) = 1550 OF Ou a= tS Se LOLG ee bigs 

19 949 -12.4 S64) 6215.49) 1010) HI 287 1045. 1158 

20 200n ea 12.5 893) 154 O45 fF 12 Fe HOTS: | eat 

21 LOD 12/5 922) = Loe Ee LO7 9M See Ow LOA ie 

ZZ 1043s -12.2 O52) | =O pe wa Ge Se 3c 35) 6 

PAS) LO7 See 12), 0 OSU aS 1. peo ee Mart Ooo Mul 2s.8 

24 BLOG weet 5 LOL ~ —l6iT> alZ0l 2S acl BihoG?” 3-130 

25 iS meets Oo) OAL =o 54 Gel ode aaa Cab 22 Jae Leta 

26 LOS eed OF OVO) L677 S26 25 A 7 

27 DES Greet. L096. LO. 9 oe 295 5 ee 15 ol 

28 T7202 26" 1638) 32358 e538 

IAS, LZ Ovo L154 = 17, Og eI od yl). 

30 L200 Rem, 9) EBD” | ta PESOS) Oe 156 

OL 13Seueee— 1340 1 12100 | M1 7s 14232 2= 26.0 

32 TSO meals. 6238 =—1) 29 4d462))  —1655 

33 Te tijiemm 3.9 2 1266  <L8sl, 500) -16.6 

34 W463 ea 3) 1294 -=18. Gee 5 39:7 -=16E5 

35 POS Mra 46: 1522) = Lome N1b7 ta Rl. 7 

36 LOO ol! 24 3505 nl OS 2 LOl6us 16.9 

yi LOOZmeaeto.G 1378) \-—1956'° 1654) ie h7 52 

38 H6GS Hevea to.8 | 1406 ~-19.9 1693 ale. 

39 PCGomee 16,2 (1434) =—2053 8 1/73 hah .8 


40 Eee 16,5: _1462>- =2027 3 -1769 1) —18..0 
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Table E-8. March Data — Fort McMurray cont'd. 
ee ee Sees Remmi A Ooi gl Mwai Je i ho 


Time (LST) 0930 1100 1300 1500 
Level Height Temp Height Temp Height Temp Height Temp 
(M) (C) (M1) (C) (M) (C) (M) (C) 
41 We9O) <= 2/0 
42 LS cain 2s 
43 1546 -21.4 
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Table.E=9.. March Data u-,.Fort.McMurray 
A aE EE a IT ST 2 a eet 


DATE: March 15/76 Location: GCOS 
Time (LST) 0930 1100 1300 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
amen A tiniest ei TonRe—.__ SAECD w *Ale B e Ch ce Mie aR a ee, 
0 320 -15.0 320 =i 2 322 — 7.0 
iW 360 =u ee} 353 =O 365 - 7.6 
Z 400 -15.6 387 rl he 409 = 21:0 
3 440 -14.1 421 =11.5 452 =m Ose 2 
4 480 -12.8 454 —l1.7 496 - 8.0 
5 pyle) -12.6 488 -11.6 539 - 8.0 
6 ees) -12.7 524; -11.3 582 - 8.2 
7 580 -12.8 5a -10.6 626 - 8.6 
8 602 aot WE) 589 -10.4 669 = 1853 
9 633 -12.4 622 -10.3 vias} - 8.5 
10 663 -12.3 750 - 8.6 
Leh 694 -12.1 800 - 8.7 
7 725 ay 843 - 8.8 
13 15) -10.0 887 - 8.9 
14 786 -10.2 930 ID 2 
15 817 -10.2 ». 4G = 9.3 
16 848 -10.3 LO, - 9.5 
17 876 -10.4 1061 =19.0/ 
18 904 -10.6 1104 = LO: 
19 923 -10.7 1148 -10.4 
20 941 -10.9 POL -10.7 
21 990 -11.2 1235 -10.9 
1g 1040 -11.4 L276 -11.0 
758) 1062 -11.7 1322 on Gt LE 
24 1085 -12.0 1365 -11.4 
25 1129 -12.3 1409 -11.5 
26 73 -12.4 1452 eltlhay/ 
27 #205 -12.7 1496 =11A9 
28 1236 -12.8 1539 ~12.0 
29 1268 ctl Ce ees 1583 -12.3 
30 1300 -13.4 1626 -12.6 
31 L331 -13.7 1670 -12.9 
32 1363 = i359 1/13 -13.1 
33 1395 -14.2 Pi5y -13.4 
34 1426 -14.4 1800 -13.5 
35 1458 -14.7 1844 -13.7 
36 1490 -14.8 1887 -14.0 
37 Les PA -15.0 1930 -14.3 
38 1553 -15.4 1974 -14.7 
Bie iksy oir -15.6 2017 -14.9 


40 1620 -15.8 2061 -14.8 


106 


Table E-9,. March Data - Fort McMurray cont'd. 


SSMS 


Time (LST) 
Level Height 
(M) 
41 1654 
42 1688 
43 ef 22 
44 igs) 
45 L789 
46 1823 
47 1856 
48 1890 
49 1924 
50 1958 
ay 1397 
a2 2025 
a2 2059 
54 ZO 
De) 2102.0 
56 2160 
By 2194 
58 2221 
oe) 2261 
60 2295 


0930 


Temp 
(C) 


-16. 
-16. 


Height 


(M) 


1100 


Temp 
(C) 


Height 


(M) 


2104 
2148 
291 
2239 
2278 
2322 
2365 
2409 
2452 
2496 
2539 
2583 
2626 
2670 
24413 
2US7 
2800 
2844 
2887 
2931 


1300 


ce 


Temp 


(C) 


-14,. 
== ere 
-15. 
malts) 
callie 
-16. 
=16). 
at Gr 
-16. 
=e 
-17. 
-17. 
—17. 
alei. 
-18. 
-18, 
-18. 
=13. 
ai), 
=A, 
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Table E-10.; March Data — Fort. McMurray, 


DATE: March 16/76 Location: GCOS 
Mie TSA Ce a ee ES 2 SASS RT 
Time (LST) 0930 1100 1300 
ee ee ee ee eee ee 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
a Se ee ee 
0 ZA = 1.8 318 - 4,6 318 0.0 
1 Sak m= O S{oulk - 4.8 358 = 0.2 
Z 421 =f) 403 ee Vas SEM) - 0.1 
3 458 - 7.7 445 - 4.0 433 roa Oieol. 
4 494 - 6.8 487 met, 469 = Ook 
5 524 ~ 4.8 530 - 3.8 499 0.0 
6 553 = 4,3 DID - 3.4 yas) 0.0 
7 589 - 4,2 614 =- 3.5 565 = O71 
8 626 = 4,3 656 - 3.6 601 - 0.4 
9 662 - 4,6 699 on Se, 637 - 0.6 
10 699 - 4.8 741 - 3.4 673 oe el) 
tal pays) - 4.6 783 = 73.3 709 oot ok ees 
12 12 - 4,5 825 - 3.4 745 - 1.8 
368: 805 - 4,6 866 = 350 782 - 2.0 
14 838 - 4,8 907 - 4.0 819 td ge 
15 871 - 4.7 949 = 4.2 855 - 2.6 
16 905 a4 990 - 4,4 So72 ~ 2.9 
5 Er 938 - 5.4 1028 "= 4,6 930 a Shy 
18 974 ee 1066 - 4.8 967 - 3.4 
19 1004 med 1105 - 5.0 1005 = 50 
20 1037 - 5.8 1143 = 5.1 1043 ~ 3.8 
Pal 1071 - 6.0 L131 - 5.4 1080 25.9 
22 1104 - 6.4 1220 = De: / 1118 elo ie: 
23 1137 = On) 256 - 5.8 1156 - 3.8 
24 1170 - 6.7 1296 ol Sipe! E1958 - 4.1 
25 1204 - 7,0 1334 = 16.50 1230 - 4.1 
26 2 - 7.2 1373 - 6.1 1268 - 4.4 
27 H275 - 7.4 1411 - 6.4 1305 - 4.9 
28 1313 a6 1449 - 6.6 1343) - 4.9 
29 (hehe e756 1488 - 6.8 1382 - 5.2 
30 590 - 7.8 1526 - 7.0 1421 - 5.4 
31 1428 - 8.0 1564 serch D 1446 ~ 5.6 
Si/4 1466 -~ 8.3 1603 mo fegt 1470 aie LS 
33 1504 - 8.6 1641 - 7.4 1499 - 6.0 
34 1543 - §.9 1679 - 7.6 1528 - 6.1 
SiS) E58 - 9.0 LIL? - 7.9 pew! - 6.3 
36 1619 - 9,2 1756 - 8.0 1586 - 6.3 
37 1658 - 9.4 1794 = Gee 1615 - 6.6 
38 1697 = D7, 1832 = Oe 1644 = Oey 
39 1736 — £99 1871 - 8.8 1672 - 7.0 
40 E77 -~10,2 1909 a GO oe) 1699 oot ry”4 


108 


’ 


\ ee a “ 
aaa 
4 ) ; bi a 
“— 


Table E-10. March Data - Fort McMurray cont'd. 


Time (LST) 0930 1100 1300 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
4l 1814 -10,4 1947 - 9,4 dy Ped | - 7.4 
42 1353 -10.5 1986 - 9.6 W755 - 7.4 
43 13915 =LO.7 2019 - 9.8 1784 - 7.8 
44 1932 -10.8 2053 -10.1 1814 = 9 
45 Loy Lk -10.8 2087 -10.3 1844 - 8.2 
46 2010 -10.8 2124 -10.6 1873 - 8.3 
47 2049 -11.0 21D -10.8 1903 - 8.4 
48 2088 =11.2 2188 -10.9 1938 - 8.7 
49 24127 -11.4 2222 ~Lied 1966 = Oeil 
50 2166 -11.8 2256 -11.3 2000 = 9,2 
yh 2205 =12.1 2290 -11.5 2033 - 9.4 
DZ 2244 = 1222 2324 -11.6 2067 - 9.6 
Do 2283 -12.5 2357 -11.8 2101 - 9.8 
54 23:25 -12.8 2598 —~12.0 2134 = 9ht9 
55 2362 -13.0 2425 -12.3 2168 -10.0 
56 2401 —13,2 2459 -12.6 2202 -10.2 
57 2440 -13.4 2497 -12.7 2235 -10.3 
58 2479 =13%°/ 2535 =13%3 2269 -10.4 
59 2518 -13.8 2574 -13.6 2302 -10.7 
60 255) -14.0 2612 -13.9 2334 -11.0 
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Table E-1l. March Data — Fort McMurray 
coi es Bea Se a a neal eR SLL. ti id ee eee a 


DATE: March 17/76 Location: GCOS 
Time (LST) 0930 1100 BESS 
Level Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) 
0 320 0.0 320 5 318 amy 
i 355 - 0.6 305 1S 355 Sys) 
Z 389 - 0.8 387 38) 393 4.9 
3 423 al LAU) 421 1.4 425 4.6 
4 457 ately: 454 . H Dee 458 4.2 
S 491 els 488 1 eye 495 4.0 
6 D2D el eS 522 1.0 ebyenh 3.6 
u 559 - 0.9 BS)! 0.8 568 3.4 
8 Bos rai 589 dle’, 605 B kegs 
9 628 0.9 622 Lae, 638 3.0 
10 662 0.9 656 TL 671 Zid: 
11 696 0.8 690 ve 705 2.4 
12 730 0.4 T23 AY 738 2s 
13 764 0.6 ViSW) eZ 770 Jib, 
14 798 0.8 790 i Bai 803 Lo7 
1S) 832 Org 823 0 835 iby: 
16 866 = O.1 857 0.8 868 2 
17 900 ai) 2 890 0.8 900 J ee 
18 934° 4H. Q52 923 0.6 335 0.9 
ee) 968 = 0,1 Seyi! 0.6 965 0-9 
20 1002 0.0 990 0.8 997), 1.6 
21 1038 0.0 1024 il 1032 20 
Paps 1074 a 24 1057 sD 1066 Zee 
72S) 1110 0.2 1093 1.6 1100 2.4 
24 1146 0.1 1130 is) 1134 Pah 
2D 1182 - 0.1 ELOY 4 1169 3.6 
26 1218 0.0 1205 ins) 1204 Beg 
rast -.+ -4254 cat ese 8 1243 Bee 1240 4.0 
28 1290 0.1. 1262 neo 1276 ore!) 
29 1332 On 1325 3 150 30 
30 13/5 O.2 13790 hs) 1346 36 
31 1417 0.0 1410 eave 1386 3.4 
32 1460 = 0.2 145105 a 1426 3.0 
39 1504 - 0.3 1491 Lies 1466 TASS 
34 1549 - 0.2 1531 a0) 1506 72a 6) 
oy 1594 - 0.2 1584 0.8 1546 IRS 
36 1638 - 0.4 1637 Ore, 1585 2.05; 
7. 1683 - 0.4 1669 0.6 1631 a 
38 728 - 0.4 1700 O25 1677 lees 
39 1773 - 0.6 ioe Oe L722 ee 
40 18h7 - 0.7 1764 On 1768 0.8 
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Table E-11, March Data = Fort McMurray cont'd. 
Se ee ee ee 


Time (LST) 0930 1100 1315 
ere 
Level Height Temp Height Temp Height Temp 

(M) (C) (M) (C) (M) (C) 
41 1868 sa OS) NAS - 0.2 1814 0.6 
42 1919 male? 1830 - 0.6 1860 0.3 
43 1970 - 1.4 1863 = 0.8 1911 Ox2 
44 2021 a L.7 1897 m= erLeer2 L963 0.0 
45 2055 = 1.9 1930 - 1.5 2021 ent AU ”4 
46 2089 eae sk 1963 nO 20:7 9E% - 0.4 
47 2123 =A ASS) 1996 - 1.4 2137 - 0.5 
43 2157 tie rd 2029 - 1.6 SIE - 0.8 
49 2191 seme 2067 - 1.6 2258 = 1.0 
50 2225 - 1.8 2106 - 1.8 2322 rod Bb 
51 2259 - 1.8 2163 - 2.0 2386 mie 
BZ 2294 meek 2221 mt 2 2449 mobs 
Do 2328 = 2.2 2262 - 2.4 FLASH) ret Bae 
54 2362 eae) 2304 eet U1, Zod 7 -~ 1.4 
oD 2396 meee 2 gil, 2388 = UZ 30 2641 ee A) 
56 2430 = Te mail 2472 - 3.1 2704 et, Abad 
D7 2464 CAE aro 2579 =e 2762 = 1.9 
58 2498 ot ES) 2687 =i}, 2820 =i Died: 
Be) 2532 - 3.7 2726 - 3.7 2865 = 2.1 
60 2566 - 4,1 2764 ae oe 2910 - 2.2 
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Table E-12. 


August Data - Fort McMurray 


SASS 


DATE: August 14/76 Location: Mildred Lake 
Time (LST) 0605 0936 1304 1435 1802 
Height Temp Temp Temp Temp Temp 

(M) (C) (C) (C) (C) (C) 
0 E53 20.0 ZO 27.6 203 
43 Noe) 1S ely 24.9 Ziel 25.8 
85 i760 13.6 PA) 26.8 PAS eg 
Ahan ey) 16.2 2 ose 26.6 24.9 
lige) HB YeS) LT .s Pie eS) 26.3 24.5 
Pals: Gee 17.6 2355 IRS) eee 
255 6.0 Ihe) Zones 25.0 2509 
297, 16.5 Ny Mag 23.0 25.1 23.6 
340 16.4 ALG Agee} 22.6 Dad eh 23.4 
382 16.4 Leo Wo 24.4 2320 
425 16.4 17.6 Zen PhS es) 22, 
467 L653 L726 2A. 2330 22.4 
510 16.3 Lowel ab ops} 7a ei 22.4 
Doe 1653 eae 21.0 22a PANES} 
595 L673 13 70 20.5 Z2e4 21.3 
638 1G. 2 Ba /ess) 20.0 Das M 21 
680 eyes) 170s. 3) LOe6 21.8 20.6 
722 ee, Lee ade HAS pase Zao 20.3 
765 Tigo) io 74 oe Zis6 20 e 1: 
807 LES 15.8 Los PA as) 1S ets: 
850 Wes 15.4 Lo =6 PEAS 19.4 
892 L652 1534 19.4 203.6 Pore 
S85 Meares: 155 ToS 20.9 1659 
O77. 16.3 15.2. LOL Z2ar20 Lows, 
1020 Geez 199 LO 9 2059 dike yaid 
1063 16.2 De, 18.8 PAVE) oo 
1105 POe2 LS 18.6 20.6 L736 
1147 LEB, Pe LS32 20533 L/S) 
1190 152.6 a Ba y// Lond WSS) 16.8 
L232 15.5 USya (8) 18.0 1953 16.6 
TPS 1B LS ek L1.9 18.9 LOS 
Ley AS) 1579 17.6 Los Teo 
1360 TESTE) Leo 17.4 LOL 15.0 
1402 15%: ies 17.4 ie Ib A4 
1445 jaye LD. 1736 17.2 14.9 
1487 AS) 15.0 17.4 16.7 14.7 
1530 14.4 14.8 Be) LG. 2 14.4 
BST 2 14.0 14.6 16.4 1.6 14.0 
1615 14.0 La 5 L632 15.3 IESE! 
Los 7 14.0 eo 14.9 LSet 
1700 13.4 137 Peo 14.6 13 oe 
1742 L305 14.3 eo 
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Table E-12. August Data - Fort McMurray cont'd. 


Time(LST) 0603 0936 1304 1435 1802 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (Cc) 
ese rere Ss hp As ois EEL OE, ood, ici ct ce 
1785 27 129 550 13.9 12.6 

1827 Bierae yt: 
1870 25k Wen 14.2 13.4 aw 
1912 a RAY 
1955 4 1 a es: Les) aw) 118 
LOO di 0) 
2040 L058 LOes Ves AER) 14 
2082 L052 
AVI 1OF2 cpl! 11.4 1.6 10.8 
2167 Se) 
2210 9.4 8.8 1036 10.7 LOS 
2292 8.4 
2295 eh) So 1Os2 LOKe2 O60 
2337: Leo 
2380 SoZ Fe oho 9.6 8.8 
2422 120 
2465 7.4 Ag) 8.6 9.4 8.0 
2507 6.5 
2550 6.9 OZ 7.6 8.3 oak 
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Table E-13. August Data - Fort McMurray 


DATE: August 15/76 Location: Mildred Lake 
Time (LST) 0605 0938 1301 1800 
Height Temp Temp Temp Temp 
(M) (C) (C) (C) (C) 
0 L7..0 2343 Zilina Zou 
43 16.5 2229 26.6 29022 
85 915.5 2223 262) 28.7 
127 LOSS 216 PAS Fe: 23.3 
170 uo 4 204 24.8 pee 
23 1 eS Zia 5 24.3 2163 
IE Ya, LS 2206 VBR 27.0 
297 1958 23.4 Zoe 2607 
340 204. 23.6 2329 265.2 
382 20.4 23.5 22.8 230 
425 203 Zoe 2205 25.3 
467 ZOGZ 250 22,0 Dip yes t 
510 19.9 22.9 21.9 24.9 
a2 1) 22s 2Uah 245) 
DDS 20.2 223 ZZ Zaa2 
638 20.6 2250 20.9 2326S 
680 20.4 PAS a | 2035 2555 
isa ZO v1 2404 20.4 Zon 
765 8) Ks) 2A 20.4 DE Ree 
807 19.4 2058 20.4 2253 
850 Oro 20.4 2051 21.8 
692 LO. 20.4 TS} e) JAN S| 
935 18.4 20m er) 2029 
o77 18,6 2062 Loe 2055 
1020 18.9 BORD 1336 PAUAY 
1063 iese rs) Oz 1335 19,6 
Lh05 iver ls 18.6 LOL 19°73 
LAAT] 18.6 Loe Nya t 13.9 
1190 18.4 LUCG AT p aes 18.6 
P2332 ig. 17.6 16:6 18,9 
27D 13.0 yaw 16,1 tS.) 
LL Rid. 16.3 L556 Ly. 
1360 Tyee) 16.5 Loe Li.3 
1402 IES) Lowe 14.9 17,0 
1445 Le 16.4 14.9 16.6 
1487 IBY Se 1539 14.7 16.4 
1530 16.6 eh 14.6 L661 
L572 16.4 155.2 14.2 ieory) 
1615 1620 15.0 tae eo 
1657 ay), Laney 14.2 Las 
1700 iB ele! 14.4 13.8 14.9 
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Table E-13. 
Tine(usT) 


Height 
(M) 


1735 
VEe27 
1870 
1O12 
tie ys95) 
1997 
2040 
2082 
Ze20 
2167 
2210 
pages 
2295 
Ja B59 
2380 
2422 
2465 
250/ 
2550 


August Data - Fort McMurray cont'd. 
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0938 


Temp 
(C) 


deo 


2 


LO ao 


10,6 


1020 


9.6 


Be) 


130 


Temp 
(C) 


aS eb 
10 
13.0 
Ne yay 
Ue ek 
13:0 
iat | 
1s 
aol 
ti 


1800 


Temp 
(C) 


14.9 
HE 


15.6 


Lo 


115 


‘ 


Voglatiy! basa rab aig sil 


ri Wey 7 ae ee fe eel lien 


eran fa ee. Sad ipa 


z pence ahem ate econo = aa 
jie f f (ey %: 
7 A - ’ 
ete alone gp ememening oe te nme 
vy 
$ oy 
= r; tan " 
* . 1 
F 
a | mans) MMe 
2 iu gE 
| mete de 
; 2 aes tA Cte - met ah 
7 aa Fa { ae erie : Cun 
i c od ee ee ,Jee eo) = / ate . 
i ‘ =" ‘ 
: » ES ; 
| +s | 
< ‘ Ae 
<4 
4 3 5 = i - Sy . - 
‘ : i f. a , [ ; i Be 
© x 
4? 
! } y ‘ 


ca, 


116 


Table E-14. August Data - Fort McMurray 


DATE: August 16/76 Location: Mildred Lake 
Time (LST) 0609 UIT 1254 1427 1801 
ee al ES oe ee SR ee nd RS i 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
Es pe a bier ee 
0 Io 222 a2 Orel. 23D 239 
43 S32 rae al Dias PIERS 2 
85 18.4 2007 20,7 21.0 271.0 
V2] 185.5 20.0 IRS 2652 26.6 
170 LOD i Be es 2569 25.8 ZO 
253 1328 iGo 25.6 2D aie 25.8 
2D 19.3 L728 Diane, 2502 ZOD 
297 Oey L726 2500 24.4 Doe 
340 197 Lijnee 24.6 Zed 24.8 
Boz 1978 EO 24.4 22.8 24.4 
425 ZOicD 16.4 24.0 Deir 24.2 
467 Zhe0 16.0 2a Dalia ZB iei, 
510 2An0 P5786 2269 ZA 23.1 
52, ANS) 1 re 22.0 20.8 pA T 
595 20.6 14.8 2s 20:3 222 
638 JA RAL 14.7 2069 20.0 PEAY 
680 19.3 14.6 20.3 19.4 Aes 
722. 19e4 aay 1a deg) Jade t ZO 
765 0 14,4 The Teal 18.7 20.6 
897 B67 14.0 Se 6 LOso 20733 
850 1623 13s 17.8 I, 1933 
892 Wow 1355 Vie 17.4 LEIA) 
935 Leto | 13.4 LOZ Le JhoAk LSe9 
O77, 18.0 ISAS) 16.1 16.8 18.6 
1020 1350 14.2 Lo. 0 16.4 18.0 
1063 Bp) 14.2 LO Lesa, 17.4 
L105 L7is8 14.3 16.8 isypy 17.0 
1147 78 14.4 UA e SUSY ais) 16.8 
1190 1S) 14.3 AN ER, 22 1550 L6n3 
T2372 17.4 1a 52 Ee 2 14.6 15.9 
L275 17.4 14.4 0 bel 5.6 
1317 ee. 
1360 ey ok: 14.4 aha 1B Shey) 14.6 
1402 LGo7 
1445 LO. 14.2 15 n2 13.4 14.1 
1487 LG 2 
1530 ibe) 1356 LAR 1235 13.5 
L572 ea. 0 
1615 13 1358 La 3 120 13.0 
1657 15.0 
1700 Ne 9 D350 13.4 Sew) Te 
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Table E-14. August Data - Fort McMurray cont'd. 


Time (LST) 0609 0927 1254 1427 1801 
NL ie id AGAIN a ae ted a ps eS I ah en NE at a a 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
Ti6> 14.0 aS) 72 MPLA US 29 BI es be) 
1327 
1870 ES.2 ALAR) Ts. P28 b4 
1912 
19:55 b2.5 ely, S 11.4 12.4 Li. 
L997 
2040 . 2.9 13 ei Ly 9s gy / 
2082 : : 
2120 ie 2 10.6 diy. bf BEA fe 
2167 . 
IBALL 10.9 oS LO.9 115.0 10.7 
2252 
2295 10.4 oS 10.4 10.5 LOD 
233) 
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2465 9.0 SZ ES 9.4 Oe: 
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Table E-15. August Data - Fort McMurray 


DATE: August 17/76 Location: Mildred Lake 
Time(LST) 0610 0928 1255 1425 1806 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
0 gD Wee, HOE eS) Ao eS L/P 
43 14.0 13.4 LZ. 14.4 16.3 
85 Poe 122% Hak Ale) PSY: 15.5 
7 1334 12.0 1 eae 3.3 50 
170 13 i 1% 36 10.4 253 14.8 
2S 12.6 LVSZ 10:23 eS 14.4 
2D ee 103.9 10.0 IEG 14.0 
297, LPR LOR6 926 a. 0 D359 
340 12,0 10.0 9.0 LO 2 13520 
382 aiey Sei) 8.8 9.4 L227 
425 14 S , 8.6 92 L2e2 
467 Ps Sete! 8.3 9.0 UREAT 
510 0 8.4 Sas S25 2 
DIZ it .0 8.0 ao 8.4 10.6 
595 LV 2 Joa 8.0 8.3 LOR 
638 aM TE: 8.0 8.0 ee 
680 Lit 733 Ves) Sad OSL 
722 LORD 7.0 bet ds6 9.0 
765 10.6 EO Phebe ie 3.5 
807 103 ae) ers) UCS In 9 
850 10.0 7.0 7.0 Tek T<9 
892 Jess) 7.2 6.8 720 TEKS) 
935 SS) Wee. Geo 6.7 Jas 
977 9.1 fAo 6.4 6.6 6.6 
1020 8.8 6.9 Gaz 629 5.8 
1063 8.6 6.6 Sy) 603 Diet 
1105 8.6 (Oye) 6.0 5.8 Del 
1147 8.3 Gs 5.9 Bie) Gut 
1190 7.9 539 59 baal 6.4 
1232 7.6 eyy/ Sal Sa0 6.4 
£275 Tie Siem 5.6 4.8 6.5 
1317 Sie! 6.4 
1360 6.8 Se3 da 3.8 6.0 
1402 4.9 
1445 6.0 4.9 4.8 Bia 5.3 
1487 4.5 
1530 (sya 4.6 4.1 2.6 4.9 
1572 SA) 
1615 6.1 4.3 Siw) aA T! 4.3 
1657 SAS) 
1700 6.0 Ae 2.9 ES) Be 
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Table E-15, August Data — Fort McMurray cont'd. 

Time (LST) 0610 0928 LDS 1425 1806 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
1785 Dao Sey) 230 he 2 eras} 

£827 236 

1870 5.3 BaD: 2.4 16,0 ps 
Li js) 

1955 4.9 Bh E.2 6 2.4 
1997 Rod 

2040 4.6 2.8 aah ee 275.0) 
2082 Pe) 

225 4.0 Zo = 2 - .8 Le7 
2167 nd 
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2507 - .4 ae 
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August Data - Fort McMurray 


Table E-16. 
DATE: August 18/76 Location: Mildred Lake 
Time (LST) 0601 0928 253 1425 L759 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
0 14.4 Bad E33 18.6 17.7 
43 La 16.0 Lie Lye > bey leat 
85 Lone 15,3 15,9 16.6 16.6 
17 13.5 14.6 pies) 1601 16.0 
170 136 VEEN) 14.8 Lond 15.6 
213 34 Ruse, er. t5aD Lay Z 
255 13.0 13.3 1839 15.0 14.6 
297 12277 130 L334 14S LEZ 
340 1232 iD 1320 TA, 7 14.0 
382 Ls? ec 1276 14.2 SIT) 
425 Ae 24 TZ 0 i2e2 aera 1362 
467 LOeS8 Be DieS6 Lond £35 
510 10.3 Lb. dee 13.6 1232 
baz on Ti4 10.5 13.4 su Ry) 
DID oD LTO es) 13.3 LD 
638 oie 1659 9.3 13.0 lee 
680 oe Lice Ok 1256 10.8 
722 Did. 11.6 D5 23 10.5 
73 aT LES Ne 123.0 10.2 
807 LOS 10.8 OG 11.6 10.1 
850 20.9 10:33 elt dee EY) 
892 0 G79 829 10.9 956 
hope 2535) Ose on 80 1005 9.4 
oe! 10.3 oi 822 10%2 5.0 
1020 LOGZ 8.7 Thege| 359 8.9 
1063 10.0 Sez eo oo 8.8 
1105 9.9 Leo Tal 9:6 855 
1147 9.6 Tao 6.9 oa Si3 
L190 Died 8.)D LO 8.6 8.0 
252 8.8 9.0 7.0 8.9 raA 
1275 8.6 ok 6.9 9 E58) 
ARS BFE 6.9 S23 
1360 eee) 8.6 6.6 Sa2 6.8 
1402 
1445 7.4 Sy 6.0 8.1 6.0 
1487 
1530 6.9 7.6 RS, eyo 6.0 
1572 
1615 Oe Tae 4.7 Tiel. 539 
1657 
1700 Dake) 6.8 4.8 6.8 5.6 
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Table E-16. August Data - Fort McMurray cont'd. 
Time(LST) 0601 0928 Naps) 1425 Ni D9 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
1785 4.6 Dag 4.6 6.7 A 
1827 
1870 4.1 Deo 4.0 Dag 4.5 
TO UZ 
1955 ate 4.4 4,0 260 3) Ree) 
LOD) 
2040 Lut S10, 35,6 4.2 ae 
2082 
2025 Deo Jaz Zio Su5 2.6 
2167 
2240 ded JAS ee: Ded Eg 
2257 26.6 
2295 4 Led bed Dia 
239)] Ze 
2380 co 0. 1.4 .8 20 Poe 
2422 
2465 - .8 8 inte L 1.6 2.0 
2507 
Zo90 -1.2 » - 4 ps5 1.8 
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Table E-17. August Data - Fort McMurray 


DATE: August 19/76 Location: Mildred Lake 

Time(LST) 0608 0931 L256 1420 1811 
Height Temp Temp Temp Temp Temp 

(M) (C) (C) (C) (C) (C) 
0 1k 37 L542 20.0 20.5 217 
43 Ibe I 14.3 Bao 1329 ZL 52 
85 Les 1366 1330 Lea y: 2045 
£27 E222 WieZ Le ay feel 20,.0 
170 1235 L235 1637 1659 LES: 
213 IAD 1230 16.4 foe 18.8 
25D 14.4 Lies 16.0 Bboy | 18.1 
297 14.5 Lrs0 1538 Nee 7 17.4 
340 14.6 TLes 155 Tat: 16.6 
382 14.4 Peg LES ARES 15.8 
425 1453 12.0 14.9 Ze] 14.8 
467 LG 2 A yal) 14:7 P22 14.2 
510 14.0 Tey 14.5 dba y/ 13.4 
SDZ 1356 Ves 14.2 1O%.9 12S 
525 ts Se Tit 9 10.1 252 
638 Shell LOS, 13.05 one PAS6 
680 rash 10.4 13.4 8.7 he ys 
722 2) 10,0 ss al Te6 10.6 
765 LAA 726 1258 dae 10.0 
807 PAG 9 35 12.4 620 9.4 
850 1 es) 9.0 2.0 By 851) 
892 eles 2 8./ 1 Es, Diet 129 
935 Oe 8.4 Lt. 3 5.0 7.4 
Oe), 1036 Sr 10.8 4.4 6.6 
1020 10.3 Leo 10.4 3.60 Oa 
1063 ee) Teo 10.0 Sa2 ie 
1105 Sis, Te! 9.4 Sha ak 4.3 
1147 oes 6.8 9.0 Jase) 3.8 
1190 Se 6.6 Say 2k Bs 
£232 Sa2 Sy) 8.4 page| Vig 
1275 fad 6.1 Sez aps: 2.0 
AES WEY Del ibe) bo 
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1402 yr 
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1615 58 Sys Det 6 Pol 
1657 £0 
1700 Diz 3.0 4.8 3 oo 
1742 Zao 6 
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Table E- 17. August Data - Fort McMurray cont'd. 


Time(LST) 0608 0931 26 1420 1811 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
1785 4,3 ZV 43 gl 4 
1827 201 we 
1870 326 1.8 3.8 ro gal = 14 
19 1.6 -'.7 
£955 DE) | Be a Beg =1.8 Lee 
1997 ined b et yd 
2040 Ne) 10 Za3 ava 
2082 16 -2.0 ev hel 
2025 a we 255 = 2e0 mee 
2167 a itie: -3.4 =2.0 
2210 26 =D Des PA 
222 can A =—3.9 =332 
2295 = il =e Zaz -2./] -3.0 

255) -2.9 

2380 = 9 -1.3 230 -3.4 
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2550 -2.3 oi 9 5h) -6.1 -4.6 


1 


3 


Ay ‘ vay aay * - 
peperin nas tole lap ab mh) Selena Nive, some 6 Gin oe 


0 iS Gee 


Table E-18, August Data - Fort McMurray 
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DATE: August 20/76 Location: Mildred Lake 
Time(LST) 0613 0943 1313 1444 1805 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
0 9.4 12.6 17.8 18.5 194 
43 Se) LLG Gs DiS LO 
85 9.4 bist LOn2 1605 LGE4 
127 ONS: LO, L528 16,0 18.0 
PavAs) 11.4 LS 1574 15.4 Py. 30 
213 7 Bi Ss) LZ 14.8. ay Aral 
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297 11.8 ue ea 15,0 La 16.0 
340 Lee Lg 14.9 13.8 15.8 
382 LEA 10.7 1425 13.4 14.8 
425 ae 2 10.5 13.6 1320 14.3 
467 10.6 10.3 ASeZ L206 JESIALS 
510 LO:..6 LG. 0 134 122 122.9 
552 10.4 9.4 aS Liss 22 
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Table E-18. August Data - Fort McMurray cont'd. 


Time(LST) 0613 0943 Pes 1444 1805 
Height. Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
L785 8 45] 4.0 Ph 8) Zak 
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Table E- 20, August Data - Fort McMurray 


DATE: )fAuecust 20/76 Location: Mildred Lake 
Time(LST) 0616 0943 1316 1441 1801 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
0 DuF L565 Nits, UG 54 18.3 
43 tO 14.6 16.4 18.8 eS 
85 1002 Lael LEAs) 1330 
27 9-9 pee Yoys Loyt 72 
170 10-9 Lowe UE araglb 1750 L635 
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1 cba d245 10.6 14.9 11.6 
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638 10.6 a4 759 13.0 8.8 
680 LOm6 Leal 20 Ze 8.0 
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807 9.9 LORS 6.3 10.8 6.3 
850 9.6 10.3 Wee 5.8 
892 9.4 10.1 952 £0 x2 D260 
935 9.0 9.0 PSO 92:8 455 
977 Sei 9.4 as WH O35 339 
1020 8.5 9 2 14.6 8.6 Lae 
1063 Cae o.9 West 8.0 
1105 7.8 8.6 
1147 hes S33 SAM. 
1190 des 8.0 6.6 
4232 6.8 ie 6.1 
1275 6.6 Tak 5.4 
1317 Bue 
1360 529 625 4.9 
1402 4.5 
1445 5 6.2 4.2 
1487 4.2 
1530 4.6 5.8 3.9 
Lp i2 oye? | 
1615 Shee) S30 S5/ 
1657 3u3 
1700 Ba, 4.8 Bw 
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Table E-20. August Data - Fort McMurray cont'd. 


Time (LST) 0616 0943 136 1441 1801 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
1785 253 She] 3.0 

1827 

1870 1.38 3.4 Cine 
19a2 

1955 1.0 Died (aso) 
L997 

2040 4 2.2 ie 
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Table E-21. September Data - Fort McMurray 


DATE: September 23/76 Location: Mildred Lake 
Time(LST) 0730 1039 V5 G2 Va 27, 1807 
Height Temp Temp Temp Temp Temp 

(M) | (C) (C) (C) (C) (C) 

0 0.6 7¢4 22 3.3 LoS 
43 Ze 324 fig ies) dfs pays TS ss 
85 PRS 4.9 Tn gs0) Jel) LOO 
id) o.6 4.7 1007 1a) 19.9 
170 Sie) 4.1 1023 ING 19.3 
2A SHsKe) 336 9°38 M8 Raa 1955 
259 Se 4.4 GSS 9.9 dbs ley 
2o7 Sig 0 Sere: 1050 18.9 
340 Thicke: Thee? OFS 9.9 18.4 
382 Sig? Te) oS Sel 18.0 
425 LOn5 O20 OTD oy Bees 
467 LOs7 a) 9275 o06 aL 
510 1004 Seal Lon ORG Lo. 7 
Sys 10%2 Om 1038 LOS Len 
ye fs) 10.1 RA LON? 1G Ler 
638 LOe0 Vhs) EGrS 10.0 Le eg 
680 9.8 ee 10.1 Sa 15.6 
ep gs) aah 1050 9.4 1586 
765 9.4 Ore O60 Onl nis 922 
807 93 Sgt 9.4 8.8 14.8 
850 oe 970 oe 8.5 14.4 
892 O22. Seo 8.9 atl 14.1 
33D 8.8 8.6 eae) 7.9 1335 
o7d Ses ahi) 942 736 ey) 
1020 ia 8.0 O2 Bes L335 
1063 -airgs) Fe) Sel FES eek 
1105 S73). Pee) 8.8 7.4 13.0 
1147 Sel Fe 7.8 124.8 
1190 Fey! 6.8 7.4 6.8 12;.6 
L232 aS 6.6 Pee 6.4 12.4 
Nay es 6.9 6.6 Tispe 653 1252 
1317 6.6 625 Tp BEA 
1360 6.4 6.6 TAGS 6.3 BES Paty 
1402 6.0 6.0 Ife ah 6.2 174 
1445 Sys) bye) 621 L039 
1487 Dino ae Bers) 10.6 
1530 Sry) 58 6.5 yA 7 10.3 
1572 5.4 Diaie Oot Bars 10.1 
1615 5.3 Dat D9 Dien o.9 
1657 Del 4.5 DD) 4.6 9.8 
1700 4.8 4,4 5.4 4.4 Shs tl 
1742 4.6 4.5 Dark 4,2 | 
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September Data - Fort McMurray cont'd. 
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Table E-22. September Data - Fort McMurray 


DATE: September 24/76 Location: Mildred Lake 
Time(LST) 0725 1043 1338 1437 1811 
Height Temp Temp Temp Temp Temp 

(M) (C) (C) (C) (C) (C) 
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297 oS 10.4 9.0 Dal 16.0 
340 2) A 10.5 Oe 10.8 1556 
382 cURL 107.5 ee) 1022 15.4 
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467 10¥7 TOS 8 LOZ als) 
510 Le: LOR 2 LOy3 1035 Is 4k 
SoZ 1igZ LI 10%) LIZ 14.9 
59D CS ea LQ: 1 LL 14.5 
638 LA 2 D2 4AN 2 11°50 L453 
680 b1..3 1253 I es 1329 
P22 ede eS: ORS. iO 13.6 
765 ABE: 12702 10.26 L029 13.5 
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1700 5.4 5.4 4.5 4.0 6.4 
1742 57.0 4.0 SN 6212 
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September Data - Fort McMurray cont'd. 
(ce) 
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September Data - Fort McMurray 
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September Data - Fort McMurray cont'd. 
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Table E-24,. September Data - Fort McMurray 


DATE: September 26/76 Location: Mildred Lake 
Time (LST) 0716 1056 W355 1440 1807 
Height Temp Temp Temp Temp Temp 

(M) (C) (C) (C) (C) (C) 
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September Data - Fort McMurray cont'd. 
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Table E-25. September Data - Fort McMurray 


DATE: September 27/76 Location: Mildred Lake 
Time (LST) 0748 1057 1330 1430 1752 
HBS Se Sa A SE 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
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Table E-25. September Data - Fort McMurray cont'd. 
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Height 
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Table E-26. September Data - Fort McMurray 


DATE: September 28/76 Location: Mildred Lake 
Time (LST) 0734 1104 AS AWS 1812 
Height Temp Temp Temp Temp Temp 

(M) (C) (C) (C) (C) (C) 
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Table E-26. September Data - Fort McMurray cont'd. 


Time(LST) 0734 1104 1341 L735 1812 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
1785 7.0 Oe Eis) 116 11.4 
1827 6.8 OE 4 LZ. 11.4 
1870 6.6 8.5 eS LOe3 Nena 
EoE2 6.4 0 10.4 10.9 
195) 6.1 Tiss, 10.6 10.0 10.5 
2997, 6.1 10,4 9.4 10.4 
2040 5.9 TZ 1031 8.8 10.2 
2082 6.0 9310 8.4. HO. 
Pas 6.1 6.6 on BeZ 10.0 
2167 6.5 8.8 7.6 oo 
2210 6.6 6.0 Soy dis Ss, 
2257 6.8 8.2 Se) 9.3 
2295 6.7 5.6 6.4 a0) 
2337 6.6 8.0 (oa 8.7 
2380 6.3 De Cn 5.0 8.4 
2422 Die) 6.8 By il Saez 
2465 51) 4.8 6.4 Dee Sa2 
2507 a 67.0 Dye. 7.6 
2550 4,8 4.0 Sas: 5.0 ee 
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Table E-26. September Data - Fort McMurray 


DATE: September 29/76 Location: Mildred Lake 
Time (LST) 0754 1109 1353 1807 
Height Temp Temp Temp Temp. 
(M) (C) (C) (C) “n6G) 
0 15,6 L6ye 7, 16.0 13.7 
43 15.6 Lon 14.6 17.5 
85 Ss, 3 L536 14d 17.0 
127 ubsiaa 552 135% 16.5 
170 14.9 14.9 12.9 15.6 
213 ns 7 Las 12.3 1523 
255 14.6 14.2 Lie 5 14.1 
297 143 13.8 10.8 1A 
340 1 3s.9 T3ns5 10.0 13.0 
382 13.64 sac ae) 9.0 Zed 
425 13.4 12.6 Souk Lk. 9 
467 £330 126) ho 10.9 
510 ee) 11.6 7.0 9.6 
SOL LPF Vier 6.35 25 
595 14 10. 5 ee) 9.2 
638 1a 72 10.0 4.9 8.8 
680 1220 9.6 4.0 6.9 
722 L.8 Dade oe) Dx6 
765 11Ls8 OL s(8 4.4 
807 220 8.8 ES 3.8 
850 B19 9.7 DAT! 34:6 
892 11.8 9.8 2d: 200 
935 11.4 9.8 1.0 2.0 
977 sal: 9.6 0.3 LAs) 
1020 10.8 9.4 -0.3 ee 
1063 LOSS oat. -0.8 0.8 
1105 10.0 8.5/7 0.9 
1147 9.5 8.5. 2 
1190 9.4 8.2 -2,2 Lez 
1232 9ez ee) 0.5 
1275 eed) 1.8 O.2 
1317 8.3 dad -3.6 Ont 
1360 8.0 6.8 -4,.1 -0.3 
1402 7.7 6.4 -3.9 -0.4 
1445 hou 6.0 -4.7 -0.5 
1487 730 5 iu) -4.5 -1.0 
1530 6.6 Sao -4.8 -1.1 
1572 Gaz 4.7 -4.8 -1.2 
1615 ene. 4.4 -4.6 -1.6 
1657 Aad -1.6 
1700 ws ik Baa, -1.7 
1742 4.8 a Ye -5.1 -2.6 
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Table E-26. September Data - Fort McMurray cont'd. 


Time (LST) 0754 1109 1353 1807 
Height Temp Temp Temp Temp 
(M) (C) (C) (C) (C) 
1785 4.4 Se) -4.7 -2.9 
1827 4.0 Bo) -4.6 =3.2 
1870 3.8 Be2 -3.9 -3./ 
1912 36D 29 ~4.5 -4.0 
1955 Oa Dh ie: -4.7 -3.9 
L997 Sia) -4.8 -4.0 
2040 Shel Zao -5.3 -4.1 
2082 350 230 -5.6 -4,4 
2125 260 Aue) -6.0 -4.4 
2167 Zed LZ -6.1 -4.8 
2240 23 0 -6,2 -4.8 
2252 2a 0.7 -6.2 -5.0 
2295 2.4 0.5 -6.2 =e 
23379 aed, O52 -6.2 -5.1 
2380 0.0 -6.0 -5.0 
2422 -0.4 -6.1 -5.0 
2465 -0.8 -6.4 =) 52 
2507 =O -6.5 -5.2 
2550 -1.0 —6.5 -5.2 
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September Data - Fort McMurray 


Table E-27. 


Mildred Lake 


Location: 


September 30/76 


DATE: 


1102 1354 1445 1808 


0756 


Time(LST) 
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Table E-27. September Data - Fort McMurray cont'd. 
Time (LST) 0756 1102 1354 1445 1808 
Height Temp Temp Temp Temp Temp 
(M) (C) (C) (C) (C) (C) 
1785 0.0 0.4 -7.6 -2.6 -2.7 
13827 =0. 2 -0.1 -7./ -3.8 -3.0 
1870 -0.4 -0.2 -8.0 -4.2 -3.6 
L912 -0.9 -0.6 -8.6 -4.5 -3.6 
1955 -1.3 =9.0 =—Deo -4.2 
1997 -16 -1.4 -9.3 =-5./ ~4.8 
2040 -2.0 -1.4 -9.6 -6.4 
2082 ob eh! -10.1 ~6,/7 
2125 -2.6 —2 eet -6.9 -5.4 
2167 Od -2.5 -10.5 -7.0 -5.8 
2210 ; -3.9 -2.8 -11.0 aiiodl -6.4 
Pape -4.0 =U -11.4 =e. -6.9 
2295 -4.2 -2.8 og Ey -8.9 -7.2 
fabio) -4.8 -3.8 bls) -9.2 -7.4 
2380 -5.3 -4.1 =e 1 eos: -7.8 
2422 -5.8 -4.3 = -10.6 -8.1 
2465 -6.1 -4.7 -12.3 lier -8.4 
2007 == a -11.5 -8.5 
2550 -6.7 -5.3 ding se edeglt 
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Table E-28. November Data - Edmonton. 


DATE: November 22/74 Location: Edmonton 
Time (LST) 0715 0828 0955 1100 
Level Height. Temp Height Temp Height Temp Height Temp 
(M) (C) (M) (C) (M) (C) (M) (C) 

0 0 -15.5 0 -15.1 0 =1207 0 -10.4 

uk 40 12 a7 39 -12.6 yi -12.7 64 -10.8 

2 124 - 8.4 133 - 6.9 154 - 7.2 167 - 7.1 

3 236 - 7.2 236 - 5.4 242 ="6,'k 25, - 6.1 

4 300 - 6.1 324 - 4,5 328 - 4,9 347 - 4.6 

5 414 - 4.0 416 - 4,9 434 - 6.2 

6 502 - 4.9 501 - 5.5 521 - 6.7 

U 588 - 5.4 6024, 27>. 6.0 606 - 7.2 

8 676 - 5.4 689 - 6.8 697 - 7.6 

9 765 - 4,8 759 - 7.4 780 - 8.1 
10 854 - 6.7 847 Nh 860 - 9.0 
Bal 944 - 7.8 935 - 8.4 948 - 9.5 
12 1033 - 9.5 1020 - 8.9 1028 -10.0 
13 2122 -10.5 1109 - 8.4 da eS) -10.1 
14 212 -11.0 1199 - 8.6 Loe - 9,3 
15 #301 -12./7 1290 = 9.1 1299 - 8.8 
16 1390 -12.9 377 mi Qi? 1383 - 9,2 
iW) 1479 -11.4 1459 -10.0 1472 -10.0 
18 1569 -11.5 £542 -10.4 1575 -10.4 
19 1658 -12.4 1628 -10.6 1608 -10.5 
20 TW) -13.0 a ie -l1.1 Lj5J -10.9 
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Table E-29. November Data - Edmonton 


DATE: November 23/74 Location: Edmonton 

ue fees SEPT Ge Se Fie eee 
Time 0800 0930 1100 1230 1400 
(LST) 


Level Het Temp Het Temp Het Temp Het Temp Het Temp 
(M) (COT Os) (cy) (C) 7" 4) CC) a Ct) (C) 


en a SS ID 


0 Oa LT4 OF 1 Le O78 5..6 Op eal ears) 0 0.9 
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4 B25 = noo Sle = 4i0- 27 Gee ON sees 0.7 |) BAZ 0.9 
5 US Gawen on O97) — 2.9" SIR TRan Coed imme Od ue oro. (— Le? 
6 Solem es Omto>) 3.7 GOAN Ee) yee oZomennO we COON — 1s.0 
7 SOS Poeb ss ola (3.8 | SO7eP =. 11s RO hom cert TSO) ett hei ie 
8 673 a5. om F000 — 4247.) 653. 716) FO2 Eee te ll 204 
2) JOSE ole — 45° aI war Onn Olea On LO cue ho 6 L 
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20 172 eto) L717” 10. 5- 7200 ne wom Oo tO. 4 20625 5-1 029 
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Table E-30. November Data - Edmonton 


DATE: November 26/74 Location: Edmonton 
Time (LST) 0800 0930 1100 1400 
Level Height Temp Height Temp Height Temp Height Temp 


(M) (C) (M) (C) (M) (C) (M) (C) 


0 O> '= S44 O° 14a OF aa 1 5 0 Et 
1 699 = G2 65." a2 65 He 26 65 O&5 
Z Lee 93.5 187 . = B.8ea, 284 56—-3 8 195m 2048 
3 290 F742 309° = 4.2 B02 Tee 465 BASE PZ 
4 402 - 4.5 431 - 4.8 420 Ba) a2 4860) — 93.6 
5) 12 = 2 Doo (be. 4 D398 Bee to eZ 59bMe = +452 
6 61S 5 -— "622 674 - 6.0 6569 5.5 685 -—- 4.6 
7 (OSM 1652 796 = 6.4 114 Se 5.36 160m -—F49 
8 738.4% —%OR9 ONG. Vga. O 093 Bre O20 S46) = )oak 
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2 UES S6987 «6405 PRED. 2 wi oSee OG UN ie eae 7 
13 1241) -10.4 TASH SEO Ne 121 = 5S 
14 P3258 4 ES 2 POOR 9 T=) 7 ela L308 = 5.0 
15 1452-119 1 20. 9Gi— 7 ANS L372 UF 558 
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L7 1639)4, 12.9 VO5 Guat ia= 7 Gn, L497 = 6.8 
18 Lb) =E3 51 2074 80 Lae =e 
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20 1920 -14.2 23) gar - 88), ) L238) — 41.6 
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Table E-31. November Data - Edmonton 


DATE: November 27/74 Location: Edmonton 
Time 0800 0930 1100 1230 1400 
(LST) 


Level Hgt Temp Het Temp Het Temp Hgt Temp Hgt Temp 
(M) Ke), ct) (Cc) (@) @) 1) (cy Gh) (C) 


————E ___________________ = a 


0) Ob naoeS Q  a8.8 Ob 258 Op 2022 0 2 
1 53), 643 68 -4.9 64 -=3.3 54-1 cy 0.4 
2 Nee 9 eo, ae) 03 RGA aS BP Ae eB | 160) 047 
3 223 =0, 6 ~261 On Ti 260) @- Ora aor Oe 25 wae, 
4 300 -0.4 349 0.4 3507, —0.9 300) ~-0.2 346 -=2.0 
5 S3e J -Or4n 6436) . -0.2 480 On 9 40S 0.3) -435r" -E.E 
6 Gly ~-=O..4. (S20 0.3) 2525) A0R9 496 -O0%.  S2u.  Ehi5 
7 SH2e HO. me | 604 0:9 Glo —lnl a 5ee. 1-0. 602) IES 
8 624 -0.5 688 OL Oy FOS NEON, 670 8%. am (686)! —be8 
9 Ode Oley R73 0:7, 796 —<0.2. 76% ~-1.G “769 1.9 

10 339) LezOnse Sere 40.11 SB O50) Saar) Ka) | 854-16 

D1 Salk 205%. 955  =—0..8) | 970) 1 —0..5. '9R6H! 7-17 950) 1-260 

12 O53) GV? 91040 =<0.9 1070) =Pek’ 20B29-2.% 1044 ~=284 

13 TOSS) Sr 123) HT Se SG 2 ay. = 26 es | =3.0 

14 itn 2S $2 ok ea aro) Ne 9H 3.3- 230) 37 
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16 W282) V2 368 «3. OB ee OO! HO) Os! “AR 

17 1364, 92384 444; ) <6 Say 5, GT) Se 1492. 9 587 

18 UihG. bad) 1586 ——5. 4." enes © 5 9) 1576), SoG (58h 685 

19 152815 24.7. W694 —Geli “1G86) Ya6.6 1655) =6.5) 1669) . -=7.0 

20 iGO: a5 is WSS. a7. 0: ee Ser laa Sah 59). HIE, 
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Table E-32. November Data - Edmonton 


eS 


DATE: November 28/74 Location: Edmonton 
Time 0800 0930 1100 1230 1400 
(LST) 


Level Hgt Temp Het Temp Het Temp Het Temp  Hgt 
(M) (Cc) (M4) (GG) 00) (GC) OD) CG) Cs) 


= 
a§ 
Vg 


0 0 -10.8 QO =-9.6 OS 789 Q ele 0 0.7 
i. 64\- =7,6 o7)) -726 63) 05 1 SL S22 Py 0.6 
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Table E-34. November Data - Edmonton 
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